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Abstract 

In this paper, we mainly employed the idea of the previous paper [36] to study the sharp uniform 
estimates with 1 < p < oo for more general elliptic systems with the Neumann boundary condition on a 
bounded domain, arising in homogenization theory. Based on the skills developed by Z. Shen in [29] 
and by T. Suslina in [33l|34], we also established the convergence rates on a bounded domain 
and a Lipschitz domain, respectively. Here we found a “rough” version of the first order correctors (see 
(11.121) 1. which can unify the proof in and [33]. It allows us to skip the corresponding convergence 
results on that are the preconditions in [33][34]. Our results can be regarded as an extension of [24] 
developed by C. Kenig, F. Lin, Z. Shen, as well as of [34] investigated by T. Suslina. 


1 Introduction and main results 

M. Avellaneda and F. Lin developed the compactness methods in [Sill] to study uniform regularity es¬ 
timates for Dirichlet problems in homogenization theory in the end of 1980s. For the Neumann boundary 
value problem, it is not until [23] established by C. Kenig, F. Lin and Z. Shen in 2013 that there was no 
significant progress on this topic. Recently, a new method has been introduced in [2l[29] by S. Armstrong 
and Z. Shen to arrive at the sharp regularity estimates, uniformly down to the microscopic scale, without 
smoothness assumptions, for Dirichlet and Neumann problems in periodic or non-periodic settings. Mean¬ 
while T. Suslina [331134] derived the sharp 0(e) convergence rate in L^(U) for elliptic systems with either 
Dirichlet or Neumann boundary conditions in a domain. 

Inspired by these papers, we originally investigated some uniform regularity estimates for the elliptic 
operator with rapidly oscillating potentials that is 

= —^Ue H— 'W(x/e)ue + Xu^ = F in U, 

where W is referred to as the rapidly oscillating potential term (see [5] pp.93]). As we have shown in [36] . 
the operator is only a special case of Ce, and therefore indicates that our results are not very trivial as 
it seems to be. 

Returning to this paper, neither the well-known compactness methods nor the new developed technique 
is rigidly used. Instead we try to make full use of the previous work in [24] . On account of these results, 
we mainly establish the uniform W^’'^ estimates (1 < p < oo), as well as the convergence rates for more 
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general elliptic systems with the Nuemann boundary condition in homogenization theory. More precisely, 
we consider the following operators depending on parameter e > 0, 

Le = —div[74(x/e)V + V{x/e)\ + B{x/s)V + c{x/s) + XI 

where A > 0 is a constant, and I = (e"^) is an identity matrix. 

Let d > 3, m > 1, and 1 < i,j < d and 1 < a,/3 < m. Suppose that A = V = B = 

c = (c"^) are real measurable functions, satisfying the following conditions: 

• the uniform ellipticity condition 

for y e M", and C = (C) G M™'', where /r > 0; (1.1) 

(The summation convention for repeated indices is used throughout.) 

• the periodicity condition 

A{y + z) = A{y), V{y + z) = V{y), B{y + z) = B{y), c{y + z) = c{y), for y G and z G (1.2) 

• the boundedness condition 

max{||l/||ioo(Rd), ||.B||ioc(Kd), ||c||ioc(Kd)} < Ki, where ki > 0; (1.3) 

• the regularity condition 

max{||^||c-o,r(Kd), ||y||co,-(Rd)} < K 2 , where r G (0,1) and «:2 > 0. (1.4) 

Set K = max{Ki,«; 2 }, and we say A G A(/i,r, k) if A = A{y) satishes the conditions (|l.ip . (|1.2p and (|1.4p . 
Throughout this paper, we always assume that 12 is a bounded domain with rj G [t, 1), and = 
—div[A(x/e)V] is the elliptic operator from |24] . unless otherwise stated. 

We focus on the Neumann boundary value problem: 

Ce{ue) = div(/) + F 
Beiue) = g-n- f 

where 

[Be{ue)Y = ni(x)l/“^(x/e)uf + ni(x)a“/(x/e)^ 

denotes the conormal derivative of with respect to Tg, and n = (ni, • • • ,nd) is the outward unit normal 
vector to dH. 

The following are the main resnlts of this paper. 

Theorem 1.1 estimates). Let 1 < p < oo. Suppose that A G YMOPM*^) satisfies (jl.ljl . (jl.2l) . and 

other coefficients satisfy (11.31) . Let f G LP{Q; F G L'^{Ll] M™) and g G B~^/P'P{dQ; M™'), where q = 

if P > 2 and q = pifl<p<2. Assume that G 1T^’^(12; ffi™') Ci is the weak solution to (11.51) . 

Then we have the uniform estimate 

ll'“e|lwi’P(Q) < C'{||/IIlp( 0) + 11^1^9(0) + lbllB-i/p.p(af2)}) (1-6) 

where C depends only on p,u;, n, X,m,d,p and Ll. 


in 12, 
on 512, 




(1.5) 
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Here we refer the reader to [30l pp.2283] for the class of VMO(M'^), and M”*) denotes the 

Besov space of order a (see IB). Compared to the proof of [241 Theorem 1.1], (where the solution belongs 
to the space of the equivalent classes in with respect to the relation — M”*), 

we can not employ Poincare inequality freely any longer. In other words, here we lack the equivalence 
between ||u£||vi/i.p(q) and ||Vue||iP(f 7 ). Thus we have to hrst estimate ||u£||vi/i’P(fJ) for p > 2 by interpolation 
inequalities, and then we can infer the quantity of ||Vtte||iP and HuglliP for 1 < p < 2 by duality, respectively. 
We remark that there are no periodicity or regularity assumptions on the coefficients of the lower order 
terms, and the estimate (II.6|) still holds when H is a bounded domain (see [29]). In recent years, 
the uniform estimates for different type of operators in homogenization theory have been studied 

extensively (see |2l[3l[Hl[l3l[Ill[2ni[2lll3ll|36] and their references). 

Theorem 1.2 (Lipschitz estimates). Suppose that A G V satisfies (|1.2I) and (|1.4p . and other 

coefficients satisfy (II.3p . Let p > d and 0 < cr < 1. Then, for any F £ LP(H;M™), / £ C'^’°'(H; M™'^) and 
g £ C°’'"(9n;]R”'), the weak solution Ue to ([I3|) satisfies the uniform estimate 


< C'{ll/llco. 


(1.7) 


where C depends only on g,T, k, X,p,d,m,r], and 11. 


We point out that (ll.7|l can not be improved even with data and domain. In virtue of the compactness 
methods, the estimate (|1.7p was established in |24l Theorem 1.2] for with the Neumann boundary 

condition under the additional symmetry condition A* = A, (that is = Oj“), while this symmetry 

condition was removed recently in |2|. By means of |24l Theorem 1.2], the proof of this theorem can be 
completed by the method analogous to that used in [36]. We hrst construct the Neumann boundary corrector 
associated with V via 


Lsi^efl) = div( 14 ) in H, 

= n - (V — Vfi) on (911, 


( 1 . 8 ) 


where djdv^ = n ■ AffijefSI, and V is dehned in (12.31) . Then yields the transformation: 
such that Ue and solve the following equations 


(Di) 


Ce{Ue) = div(/) + F 
Se{Ue) = g-n - f 


in H, 
on cin. 


and 


(D2) 


Le(Ue) = div(/) + F 


dVe 

dVe 


= g-n-f 


in H, 
on dLl, 


respectively. Then the result of [Ml Theorem 1.2] can be directly applied to (D 2 ). Due to the same reason 
as explained in [36], we need to derive 


||^£,o - < Celn(ro/e + 2), ro = diam(ll), 

which follows from the decay estimates of Neumann matrixes dehned in |24l pp.916] (see Theorem 14.2p . as 
well as 

which are the main conclusions of Corollary 14.61 and Lemma 14.81 where 0 < (Ti < (T 2 < 1 are independent of 
e. The above two estimates together with (II.6p guarantee that the right hand side of (D 2 ) can be uniformly 
bounded by the given data in Theorem 11.21 

Note that the right hand side of (D 2 ) which involves div(/) is, as a matter of fact, more general than 
that in [241 Theorem 1.2]. We hnd a simple way inspired by [31] to derive the uniform Lipschitz estimate 
for the weak solution to Lfiufi) = div(/) in D and due/du^ = —n ■ f on dLl. The key ingredient is to 
construct the auxiliary functions {vi;,k}k=i^ which satisfy 

Leive,k) = 0 in D, dve^k/due =-nkl on dLl 
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for k = 1, • • • ,d, where is the fc’th component of the outward unit normal to dQ. This is of help to 
eliminating the bad term S/xN^ix, z)nk{z)dS{z) in (|4.3I) (see the proof of Lemma ITT]) , where denotes 
the Neumann matrixes associated with Lg. Then the desired result: ||< C'||/||c 0 ’'^(r 2 ) follows 
from |24l Theorem 1.2] at once. We remark that this argument also works for common elliptic operators, 
provided that there are some decay estimates of corresponding Neumann matrixes previously. At the end of 
the paragraph, we mention that both the compactness method in [3l[31[2l] and the new argument developed 
in [2(129] may be still valid for the Neumann problem (II.5p . however it will be quite complicated as compared 
with our method. For more references on this topic, see [2V[41 fT^[20l[Ml[25( [ ^l311136] . 


Theorem 1.3 convergence rates). Suppose that the coefficients of satisfies (II.ip . (11.21) and dop. 
Assume that Ue,uo are the weak solutions to 


(He) 


J^eiUs) = F 


in n, 
on cin. 


(Ho) 


•Co(^io) — F 
Ho(uo) = g 


in H, 
on dit, 


^ Fe{Ue) = g 

respectively. We obtain the following two results: 

(i) if Q is a bounded domain, and F G and g G M™) are given, then 

ll^^e “ '“o||L2(r2) < Ce{\\F\\x,2f^Q^ + 11511^1/2,2 (an )}; 


(1.9) 


( 1 . 10 ) 


(ii) 


let Q be a bounded Lipschitz domain. Suppose that F G L'^+i(n;M™) and g G and A 

additionally satisfies A* = A, then we have 


Ue - UoWl^in) < Cs2{\\F\\^^^^^ + ||5||l2(9C!)}! 


( 1 . 11 ) 


Moreover, if uq G then \\ue — UQ 


2d 

{n) 


< C'e||uo||// 2 (f 2 ), where C depends only on g,K,m,d and 11. 


We note that M”^) = M™'). Here Cq and Bq from the homogenized equation (Hq) 

are defined in ()2.9jl . As mentioned before, we find a new type of the first order corrector £Xk{x/E)<pk with 
A; = 0,1,... , d, which together leads to 


wl = ul 


-uT, - e 




hi 


{xle)ffil, 


k=0 


( 1 . 12 ) 


where ip = (ifl) G and {xk}t=o correctors associated with Ce (see (12. ip and (12.2p ). Then 

we acquire the “first but rough” estimate of ||m£||_H'i(f 2 ) in Lemma 15.21 bv energy methods (where we borrow 
the idea from T. Suslina in (HU pp.3475-3477]), and the next step is therefore reduced to choose suitable ipk 
to obtain the convergence rates under the different conditions. In author’s point of view. Lemma 15.21 olavs 
a role as the bifurcation in the whole proof of Theorem 11.31 

Precisely speaking, when H is a bounded domain, following the arguments developed by T. Suslina 
[HIM], it is not hard to derive ||'U’£||_H'i(f 2 ) = 0{e'i) by choosing (p^ = 5e(Mo) and ipk = Seffi/kUo) (see Remark 
EH). But we fail to obtain = 0(e) due to the bad estimate (12.261) for the boundary terms. So, we 

turn to construct by choosing (po = SeiipAeUo) and ipk = k'^o) to avoid handling the boundary 

things (see Lemma l5.3p . then by duality argument, we can arrive at the estimate (jl.lOl) (see Lemma 15.51 we 
mention that our proof actually relies on the new development in [29] i. 

For a bounded Lipschitz domain, due to the skills improved by Z. Shen in [29|, we set ipo = S‘^{ip 2 e'ao) and 
ipk = S‘l{'4>2e^kUo) (see Lemma l5.6|) . and similarly obtain ||m£||_H'i(r 2 ) = 0(£^), which leads to the estimate 
(|l.lll) . The progress is that we additionally employ the radial maximal function coupled with non-tangential 
maximal function (see (12.371) and (12.381) 1 to analysis the boundary behavior of the solution uo to (Hq). We 
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remark that the thinking is originally arose by C. Kenig, F. Lin and Z. Shen in [22]. If u € 

using duality method again, we can derive the sharp estimate ||rce|| = 0{s). Although (11.111) is not 

sharp, it opens a new door to reach the Rellich estimate (see Remark 15.111) . We mention that under some 
conditions, it is possible to remove the assumption of uq £ by using a subtle technique, and we 

will show that in another paper. 

Here Ss,Se and '02e are defined in (|2.20p . (I2.24p and (I2.36p . respectively, uq is an extension function of 
tig on M'^. Thus to varying degrees, the two results with different preconditions between [Mj and [22] may 
be reduced by Lemma [52] to figure out some proper first order correctors. Note that both (11.101) and (11.111) 
do not require the smoothness assumptions on the coefficients of Cg, where the estimate (11.101) is sharp. 

We also remark that the counterpart of Theorem 11.31 in [2][29] performs a crucial role as the start 
point for various sorts of uniform regularity estimates, which marks the new way to regularity theory of 
homogenization problems. We refer the reader to [21[5l- [71[^[16lll7ll201123ll28ll29ll331[36] and references therein 
for more results. The assumption of d > 3 is not essential but easy to organize the paper, since we usually 
have different way to handle the corresponding problem in the case of d = 2. 

This paper is organized as follows. We introduce some definitions, symbols and remarks in Section 2. 
The proof of Theorem II. II is shown in Section 3, and the uniform Lipschitz estimate is studied in Section 4. 
Finally, we discuss the convergence rates in Section 5. 


2 Preliminaries 

Define the correctors Xk = iXk^): 0 < k < d, associated with as follows: 


for k = 0, and 



' Li{Xk) = div(F) 

in W^, 


1 


) and J Xkdy = 0 

(2.1) 


Ti(xf + A) = o 

in 




and Xkdy = 0 

(2.2) 

(0,1]'^ ^ and 

dkper(Y]RA denotes 

the closure of C“^(y;M™') in 


Note that is the subset of which collects all T-periodic vector-valued 

functions (see pp.56]). By asymptotic expansion arguments, we obtain the homogenized operator 


Co = -div(AV + V) + BV + C + XI, 


where A = (d“- ), V = {V^), B = {B^) and c = (c"^) are given by 


r 


Sx'I 


'Y 




a;y ^e^]dy. 


Bf = [Bf + Bf c‘‘l> = ^ [c“'5 + Br^]dy. 


3 

7/3 


(2.3) 


Definition 2.1. Let / = (/“) G F G L^(D;R™) and g G B~^/‘^^‘^{dCL-,RF'). We say G 

id^(D;M™') is a weak solution to (II.5p . if satisfies 


= - f /" • V(/>"dx + f F"(/)“dj:+ < gA 

J U 


(2.4) 
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for any (j) £ H^{Q-,W^), where 
Be[Ue,(p] = 


afU I T/up , 

-h K u 




■dx + 


/3 


.0/3^ 

dx 


^up _|_ _|_ X'lQ^cl/^dx (2.5) 


/n I dxj ^ i 9x 

is the bilinear form associated with C^- 

Remark 2 . 2 . Choose 0“ = 1 in (j2.4p . and then we have the compatibility condition 

{Bf{x/e)^ + c‘^^{x/e)u(^dx + \j^u^dx = j^F‘^dx+ < 5",! >s-i/2,2(at2)xBi/2,2(aQ) (2.6) 

for a = 1,... ,m, which implies the counterpart of ()2.6p in |24] since B = 0,c = 0 and A = 0 there. 
Remark 2 . 3 . We similarly define the bilinear form Bq associated with Cq as 


B, 


)[uoA] ^ 


T>a/3 j) d(f)^ 
dxi ^ * V dx, 


-dx + 


Li^‘ 


0 


0.0 9^0 

dxi 


+ + \u\ 




^ U^dx 


(2.7) 


for any uq, € 77^(0; M™). 


We remark that not only the solutions approaches to uq weakly in but also the flows 

converge, i.e. A^Vus + V^Us AVuo + Vuq, and B^Vu^ + (cg + \I)ue BVuq + (c + XI)uo weakly in 
L^(0;M™') as e —> 0 (see [211 PP-31])- In other words, for any fi e we have Bi;[ue,4>] —?> Bo[uo,</>] 

as e —>■ 0. Thus by Definition 12.11 we have 


Bo[uo,</>] = - / /"•V(/)“dx+ [ F"(/>"dx+< 5r,(/) > 5 -i/ 2 . 2 (as^)xBi/ 2.2 


(an) ■ 


( 2 . 8 ) 


Let (j) € Rd(D;M™') be an arbitrary test function in p2.8p . and integrating by parts we can derive £o(u) = 
div(/) + F in D. Then we return to p2.8p with any fi € iL^(D;M™), and obtain n- (^Vmq + = ff on dfl. 

By writing Bo(uo) = n ■ (^AVuo + Vuq), we can express the homogenized problem related to the Neumann 
problem (II. 5p as following: 

f Co{uo) = div(/) + F in D, ^ 

I Bo(no) = g on dfl. 

We refer the reader to |5l pp.l03] or [211 PP-31] for more comments. We also point out that Cq is still an 
elliptic operator and if A* = A, then we may have for any ^ = (^f) € 

(see [3 pp.23]). 

Definition 2.4. Define the adjoint operator T* as 

Cl = -div|7l*(x/e)V + R*(x/e)| + V*{x/s)V + c*{x/s) + AI, 
while the corresponding boundary operator becomes 

B* = n • [A*{x/e)V + B*{x/£)]. 

Furthermore, the related bilinear form is given by 


b: 


I”- = L {“SS +""+ L {F 


al3^fi£_ „a0 a , 
^ dXi ^ ^ ^ 


for any U£,(/> € R^(D;M™'). 












Remark 2.5. If G H^{Q;W^) are two weak solutions to 
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Ce^Us) = div(/) + F in n, 
Be{ue) = g-n-f on 511, 


C*{ve) = div(/i) + H in n, 

B*{ve) = b — n ■ h on 5n, 


respectively, then we have the second Green’s formula 

< > < U^i F^i^Vg:') >= < > + < u^, >, 


( 2 . 10 ) 


by noting that B£[ue, = B*[ue, Ug]. Moreover, if /j, F,H^ L?‘{Fl-, M”*) with i = 1, • • • , d, and h, b, g vanish, 
then we have 

/ u^Hdx = — / / • Vugdx + / Fvsdx. (2-11) 

Jq Jq 


Remark 2.6. To handle the convergence rates, we define some auxiliary functions via 


d 


d 


i>s(y) =«- “37(9) - ‘hf(v%ixn. c(9)= vr - vriy) - y?f{y)-^ixn. 


( 2 . 12 ) 


and 


d 


AC = wr ■■= c - Br(s) - C(!')^{xr} ^ 

AC = wr := e"“ - c““(9) - A{C} 


/y 


'&f{y)dy = 0 , 


in 


j^dQ^{y)dy 


(2.13) 


= 0 . 


We mention that the existence of dk is given by m Theorem 4.28] on account of j-y {y)dy = 0 for 
A: = 0,1,... , d. Furthermore it is not hard to see that d^ is periodic and belongs to 


Lemma 2.7 (Cacciopolli’s inequality). Suppose that A satisfies (jl.ip and (jl.2p . and f is a Y-periodic 
function in L‘^{Y;W^'^). Let x be a weak solution o/div(^Vx + /) = 0 m M'^. Then, for any q G K™' and 
B C 2B with 0 < r < 1, we have 


{£\Vxfdy) ^ (£ lx - gfdy)+ C(£ \ffdy) 


(2.14) 


where C depends only on g, m, d. 


Proof. The proof is standard, and we provide a proof for the sake of completeness. Let u = (y — q)<fi‘^ be a 
test function with any q G M"*, where (j) = l in B,4> = 0 outside 2B and |V(/>| < C/r. Then we have 

[ 4’'^AVxY'xdy+ [ 4>AS/xY'(j){x-q)dy= [ f ■ Vxdy + 2 [ 4>f ■ V(j){x - q)dy. (2.15) 
jRd jRd Jm<i JRd- 

By using (HI) and Holder’s inequality coupled with Young’s inequality, we derive 

T 
4 

and then have 


<(>"|Vx|^dy <C<! I \VcP\^\x-q\^dy+ [ cP^\f\^dy}, 

Rd jRd 

■ [ \x-q\'^dy+ [ \f\‘^dy\. 

J2B J2B J 


[ \Vx\^dy<C 
Jb 


This implies the estimate (|2.14l) . and we complete the proof. 


□ 
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Lemma 2.8. There exist € Hp^^{Y) with k = 0,1,... ,d, such that 

= “<i 1^1 =-E'S',, (2-16) 

where 1 <i,j < d and 1 < a,'y < m. Moreover if Xk is Holder continuous, then € L°°{Y). 

Proof. We provide a proof here for the sake of completeness. Based on the observation that 

£ Kk iy)dy = 0 and ^{^ 7 } = 0> (2-17) 

(The first equality of (I2.17j) follows from (12.,ip . and second one follows from (j2.ip and (I2.2p .l we can construct 
auxiliary periodic function satisfying 

=Kk ™ iy)dy = (2-18) 

Note that the equality ()2.17l) guarantees the existence of in i?pg^(y) (see [13 theorem 4.28]). Let 

obviously EJJ^ = Moreover we have 

by noting the fact that ^{^7} = 0 in implies = 0. From Liouvill’s theorem, it follows that 

4^{dik} = C, and therefore we have = 0. 

We now turn to prove € L°°{Y) based on Xk £ C°’°'(M'^). Due to Cacciopolli’s inequality (I2.14h . 

/ \^Xk\^dy<^[ \xk-Xkix)\^dy + Cr‘^<Cr'^^+'^~^ 

J B(x,r) X J B{x,2r) 


for any r £ (0,1) and x ^Y, and this leads to 



1 + r2-2 


} < C'r2‘^-2^ 


(2.19) 


Without loss of generality, we may assume Y is centered at 0. Let tp £ C^{2Y) be a cut-off function 
satisfying (^ = 1 in B(0,2/3), ip = 0 outside i?(0,4/5), and |V<y9| < C. We thus localize the equation (I2.18P 
as A(^ip6fj^) = ipE^^ -|- 2'S/ipV6'^j^ + AtpO^^ in M'^. In view of the Newton potential, we acquire 


||V0“'^||ioo(y) < 




\&y 


I2Y \x-y\ 

OO 




k=0 


k=0 • 


|a;-y|<2-''+i 


~^<\x—y\<2~^+^ \x y\ 


iy \ 

d-l J 


ICP*) 


1/21 


< C{1 + ^2-'^7 < OO, 
/c=0 


where the estimate |VG(3:, y)\ < C\x — y\^~'^ for Newton potential G is employed in the first inequality, and 
then we use Holder’s inequality in the third inequality and the estimate ()2.19p in last one. This completes 
the proof. □ 
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Definition 2.9. Fix € C^{B{0, 1/2)) such that C ^ 0 C = 1- Define 

SeiDix) = f * Cs{x) = [ f{x-y)Ce{y)dy, (2.20) 

where Qix) = e~'^({x/e). 

Lemma 2.10. Let f £ L^(M'^) for some 1 < p < oo. Then for any h £ Lper(K'^), 


ll^(•/£)'S'e(/)||LP(Rd) < C\\h\\ip(Y)\\f\\LP{M.d)^ (2-21) 

where C depends only on d. 

Lemma 2.11. Let f £ VF^’P(]R'^) for some 1 < p < oo. Then we have 

l|5s(/) - /IIlp(r^) < C£||V/||ip(„.), (2.22) 

and furthermore obtain 

I|5s(/)||l2(k.) < and ||S',(/) - < Ce^/^V (2.23) 

where q = and C depends only on d. 

Remark 2.12. The proof of Lemma [2.101 and 12.Ill can be found in [29l pp.8]. If we define to be 

Se{f)ix) = j- f{x- £z)dz (2.24) 


for any / £ L^(]R‘^) (which is referred to as the Steklov smoothing operator originally applied to homoge¬ 
nization problem by V.V. Zhikov in [37] and further developed in (MllMI)) then the estimates (12.211) and 
(|2.22|) are still true for this kind of operators (see (HU pp.3459]). By definition, it is clear to see that the 
operator defined in [29| plays the same role as the Steklov smoothing operator in [M], but seems more 
refined in view of the extra property ()2.23p it satisfied. 


Remark 2.13. For simplicity of presentation, if / is a periodic function, we will denote /(x/e) by fe{x). 
For example, we usually write ^^(x) = A{x/£) and Xk,e{x) = Xk{x/£), and their components follow the 
same simplified way as well. 

Lemma 2.14. Let Q be a bounded domain, and be given in (|2.24l) . Suppose that u £ 
and f £ Lpg^(Y). Then we have 


ll/e<S'e(u)||i2(f7\Se) ^ C£'i\\f\\L2^Y)\\u\\Hi(T^dy (2.25) 

and if f £ then 

||/£*5£(u) 11^2(9^7) < C£~2\\f\\^i^Y)\\u\\H^(s.d) (2.26) 

where C depends only on m,d and Ll. 


Proof. The proof of (I2.25p will not be reproduced here. We refer the reader to [341 Lemma 3.3] and its 
references. We now prove the estimate (I2.26p . Let g £ C'q(M'^;M'^) be a vector field snch that (p, n) > c > 0 
on dLl, and then 


1 


||/£‘S'e(^)|lL2(aQ) — / < 

C Jg^ 


g,n> \ fsSs{u)fdS = - [ diY{Q\feSe{u)\‘^)dx 
c Jn 

j^divQ\feSsiu)\‘^dx + £ < Q,V{feSeiu)) > f^S^{u)dx'^ 


< ||/e5£(tl)||^2(Kd) + e ^||(V/)£S'e(u)||^2(Rd)||/eS'£(n)||£,2(Kd) + ||/eS'£(Vn)||^2(Kd)||/£5e(n)||j;^2(]Rd)}, 
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where we use Cauchy’s inequality in the last inequality and {Vf)e{x) = V/(x/e). Since the estimates (j2.2ip 
and (|2.22l) are still valid for (see Remark 12.121) . we obtain 

\\feS£{u)\\j^ 2 ^QQ-^ < C*! ll/lliZ^y) +e l|V/llL 2 (Y)||/||l, 2 (y)||u||^ 2 (]jjd) + ||/|lL 2 (y)||Vu||£, 2 (-]gd)||u||£, 2 (]jjd)} 

This implies the estimate ()2.26p . The proof is completed. □ 

Remark 2.15. Throughout the paper, let B[P,r) denote the open ball centered at P of radius r, and the 
symbol tq only represents the diameter of Q. Since dQ € there exists R such that for each point 

P G on there is a new coordinate system in obtained from the standard Euclidean coordinate system 
translation and rotation so that P = (0,0) and 

B{P,R)nn = B{P,R) n {{x',Xd) G : x' G and Xd > 4>{x')], 

where (j) € is a boundary function with (/)(0) = 0 and ||V(;(!)||(^o,77(Rd) < Mq. Note that the pair of 

(ry, Mq) indicates the boundary character of 12. To describe boundary estimates, we need more notation: let 

T»(P,r) = R(P,r)nn, A(P,r) =R(P,r)nn. 

and kD{P,r) = D{P,kr) and kA{P,r) = A{P,kr) with k > 0. We usually omit the center point and the 
radius of B{P, r),D{P, r), A(P, r) without confusion. In the paper, saying a constant C depends on 12 means 
this constant involves both (r],Mo) and |12|, where |12| denotes the volume of 12. 

In the following, we introduce the Schauder estimates for “classical” Neumann problem. For this purpose, 
we set L = —div(^V) and C = —div(^V + R) + BV + c + AI, where the coefficients A, V, B, c do not 
dependent on e. Besides, du \ dv = n • A\/u denotes the conormal derivative associated with L, and 
B{u) = n ■ [AVu pVv) represents the same thing for £. For simplicity of presentation, we define 

P(r; 22; A-,F-,f]g) = r(^j- \F\^dx'^ ^ + ||/||l°°(z)) + + I| 5 'IIl°°( 2 A) + (2.27) 

where / G C''^’°^(12; with a G (0,1), F G L^(12;M™') with 1 < p < oo, and g G C'°’'’'(12;M™). Note that 
D = D{P, r) and A = A(P, r) for P G 512. 

Lemma 2.16. Let LI be a bounded domain. Suppose A satisfies (jl.ip and (11.41) . and f G C'*^’°'(12; 
with a G (0, r] and F G L^(12;M'^) with p > d. Let u be a weak solution to L{u) = div(/) + F in D{Q,^r) 
and dujdv = g — n ■ f on A((5,4r), where Q G 512. Then we have the following boundary estimates: 

(i) the Schauder estimate 

[V«]co,.p) <Cr-‘^|(^^|Vupdx)^+77(r;2Z2;2A;P;/;5)|; (2.28) 

(ii) the Lipschitz estimate 

II Vullioo(^) < c| I Vul^dx) " + 77(r; 2D- 2A; P; /; < 7 )|, (2.29) 


where P(r; 2P; 2A; P;/; p) is defined in ()2.27p . and C depends on pL,T, K,m,d,o-,p and the character of LI. 
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Proof. The estimate (12.2811 can be found in m Theorem 5.52], and we thus omit the proof. For (ii), we can 
straightforward derive (I2.29P from (|2.28ll as follows. For any x,y £ D, we have 

|Vu(x)| < \Vu{x) - Vu(y)| + \Vu{y)\ < r^pVujc-o.^p) + |Vu(y)|. 

Integrating both sides with respect to y on I? and divided by \D\, we arrive at 

\Vu{x)\ < r"[Vn]co,.(z 5 ) \D\-^ + (£ | Vn|dy) < c| | + 7e(r; 2D- 2A; F; /; y)|, 

where we use Holder’s inequality in the last inequality. This implies the desired estimate. □ 

Remark 2.17. Let u be the weak solution to L{u) = div(/) + F in H and du/dv = g -n- f on clH, and the 
corresponding assumptions are given in Lemma [2.161 Based on the estimates (|2.28p . (I2.29P and |18l Theorem 
5.14] (the corresponding interior estimate), it is not hard to derive the global estimates: 

max{||VM||Loo(Q), [VM]co,a(f 2 )} < C'{||/||co,a(f 2 ) + ||F]| 2 ,p(o) + |bllco.'^(9Q)}) (2.30) 

where C depends on g,T,m,d,p,a and H. In fact, under the assumptions of Lemma (|2.16l) . It is not hard 
to derive estimates for 2 < p < oo (by using the methods developed in [12]) that 

l|VM||LP(n) < C{||/||iP(Q) + ]|F|| + ll5llB-i/p.p(9r2)}- (2-31) 

We finally mention that the estimates actually hold for 1 < p < 2 by the duality argument. 

Lemma 2.18. Let Q be a bounded domain. Suppose that the coefficients of C satisfy (11.11) . ()1.3p and 
m, and f G with a G (0,r] and F G LP{Q;W^) with p > d. Let u be a weak solution to 

C{u) = div(/) + F in D{Q, 4r) and B{u) = g-n- f on A{Q, Ar), where Q G dQ. Then we have the following 
boundary estimates: 

(i) the Lipschitz estimate 

IIVullioop) < c| i \ufdx) ^ + F(r; 2F; 2A; F; /; y)|, (2.32) 

(ii) the Schauder estimate 

[Vn]co,-T(£,) < ^ + F(r; 2F; 2A; F;/; y)|; (2.33) 

where F(r; 2F; 2A; F;/; y) is defined in ()2.27p . and C depends on p,T, K,m,d,o-,p and the character ofPt. 

Proof. There are several methods to derive the above estimates. We plan to show the corresponding global 
estimates and then employ localization argument (introduced in [361 Remark 2.11]) to establish the desired 
estimates. We rewrite C{u) = div(/) + F as L{u) = div(/ + Vu) + F — BVu — (c + XL)u in H, and 
B{u) = g — n ■ f as du/dv = g — n ■ {f + Vu) on dPL. Then it follows from (I2.3n|) that 

||Vu||^oc(f,)<C{||/|bo. ^{n) + ll-^llLp(n) + \\9\\c0’'^{dn) + Ikllc0’‘^(f2) + lklltyi’P(!^)} 

< {\\f\\c°’^{n) + ll-^llLP(n) + Il5llc0’‘^(an) + + IbllvFi’P(n)} 

< C{\\f\\cO’'^{n) + ll^llLp(f2) + + Iklltvi’P(Q) + lbllL°°(Q)} + 2 

where we use the interpolation inequality in the second inequality, and Young’s inequality in the last one. 
This implies 

l|V'u|boo(f^) < C'{||/||c-o,a(Q) + ||F|bp(f^) + ||y||co.‘^(9C!) + Ibllwi.p(n)} 

< C'{||/||c-o,a(Q) + ||F|bp(f^) + ||y||co.‘^(9r2)} 


(2.34) 




















12 


where we use the Sobolev embedding theorem ||u||x,oo(q) < (for p > d) in the first inequality, 

and estimates (which follows from (j2.31l) by using the same idea as in the proof of Theorem 1.1, and 

we thus omit details here) in the last one. Moreover, combining (I2.30h and (I2.34j) . we have 

< C{\\f\\co,<T(^Q) + ll-^lli,p(o) + WdWcO’'^(an) + lkllvi/i.°°(o)} 

< C{\\f\\c°’<^(n) + ll-^llLp(n) + llfi'llco>'^(ac2)}- 

Finally, by the localization argument (shown in [36( Remark 2.11]), we can derive (j2.32p . (I2.33p from (|2.34p 
and (|2.35l) . respectively. The proof is complete. □ 

Remark 2.19. The following notations will be used frequently: 

• distance function 5{x) = dist(3:, (912), where x G 12; 

• average of function f e = f {x)dx = Je f{x)dx, where FI is a subset of and the subscript of 
fE is usually omitted; 

• boundary layer 12 \ S,., where = {x G 12 : dist(x, (911) > r} with r > 0; 

• cut-off function ifjr (associated with S^), satisfying 

ipr = I in T, 2 r, V'r = 0 outside Sp, and {S/ijjrl < C/r] (2.36) 


• level set 5,. = {x G 12 : dist(x,(912) = r}; 

• internal diameter roo = min{dist(x, y) : x,y G 912}, and layer constant cq = roo/10. 

Remark 2.20. For 0 < r < cq, we may assume that there exist homeomorphisms : (912 ^ 9Sp = Sr 
such that Ao(<5) = Q, l^r(<5) - ^t{P)\ ~ Ir - t\ + \Q - P\ and |Ap(Q) - At{Q)\ < C'dist(Ap(Q), S't) for 
any r > s and P,Q € 912 (which are bi-Lipschitz maps, see [22l pp.l014]). Especially, we may have 
maXj,g[o^co]{||V^r||i,°°( 9 r 2 ); ||V(A“^)||^oo(aQ)} < C{r],Mo). For a function h, we define the radial maximal 
function A4{h) on 912 as 

A4(/i)((5) = sup {|/i(Ap(Q))| : 0 < r < Co} V Q G 912. (2.37) 

We mention that the radial maximal function will play an important role in the study of convergence rates 
for Lipschitz domains (we refer the reader to |22| for the original thinking, and we also refer the reader 
to |26l Theorem 5.1] for the existence of such bi-Lipschitz maps). 

Definition 2.21. The non-tangential maximalfunction of u is defined by 


{u)*{Q) = sup{]tt(x)| : X G r7VQ(Q)} V Q G 912, 


(2.38) 


where = {x G 12 : ]x — QI ^ is the cone with vertex Q and aperture Nq, and Nq > 1 is 

sufficiently large. 

Remark 2.22. Let h G L^{Tl) with 1 < p < oo. For any r G (0,co) (cq and Ar are given in Remark 12.19]) . 
we can show the estimate of ||/i||LP(n\Sr)- By (12.371) . we note that /i(Ap(x)) < A4(/i)(x) a.e. x G 912 for all 
r G (0,Co). Then 


f \hfdx 

J S7\Er 



0 JSt=At{dn) 


\h{y)f dSt{y)dt 


0 Jan 


\h{At{z))\P\VAt\dS{z)dt <Cr [ \M{h)\PdS < Cr [ \{h)*\PdS. 

Jan Jan 


(2.39) 
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where C depends only on p and the boundary character. We note that the first equality is based on so-called 
co-area formula (I2.40p . and we use the change of variable in the second one. Besides, the first inequality 
follows from Remark 12.201 In the last one, it is not hard to see M{h){Q) < {h)*{Q) by comparing the 
Definition ()2.38p with ()2.37l) . 

We now explain the co-area formula used here. Let Z{0;r) = {x € D : 0 < 6{x) < r}, then Z(0;r) = 
D \ Tir- Here we point out |V(i(x)| = 1 a.e. x € D without the proof (see [HI pp.l42]). In view of co-area 
formula (see [m Theorem 3.13]), we have 


' S7\Sr 


\h\^dx = 


/Z(0;r) 


\h\^dx = 


0 J{x(^Q-.5(x)=t} 



dU^-^dt = 


0 JSt 


\h\PdStdt, 


(2.40) 


where 5^ = {x £ D : 5{x) = t}, dH^ ^ is the {d — l)-dimensional Hausdorff measure, and dSt = dH^ ^{St) 
denotes the surface measure of St- 


Lemma 2.23. Let LI he a Lipschitz domain, and A4 associated with cq is defined in Remark \2.20[ Then for 
any h G H^{LL), we have the following estimate 


ll-^(^)llL2(an) < C'||^||/fi(n\Sco)j (2-41) 

where C depends only on d,co and the eharacter ofLl. 

Proof. The proof is based on the fundamental theorem of calculus and definition of the radial maximal func¬ 
tion, and our proof follows the idea from |22l Proposition 8.4] (they actually proved a weighted inequality). 
For any 0 < s < t < cq and P G dLl, we first have 

h{At{P)] - h{A,{py) = l‘ f{HAr{P))}dr 

and then 

fCO 

\h{At{P))\<C \Vh{Ar{P))\dr + \h{AfiP))\. (2.42) 

40 

Integrating both sides of the above inequality with respect to s, and then divided by cq, we arrive at 

/■CO / fco >1/2 

\h{At{P))\ <C {\Vh{Ar{P))\ + \h{Ar{P))\)dr < {\Vh{Ar{P))\^ + \h{Ar{P))f)drj . 

This implies 

/■Co 

\M{h){Q)\^ < C / (lVh(A.(P))|2 ^ |/j(A,(P))|2)d^. 

Jo 

Hence we obtain 

[ \M{h){Q)\‘^dS <C [ ° [ {\S/h{AriP))\‘^+ \h{Ar{P))\^)dSdr <C [ + \u\^)dx, 

Jan Jo Jan 4c\Scq 

and this gives the desired estimate (I2.4ip . □ 

Let N^(x,z) = [Nfi'^(x,z)] denote the Neumann matrix associated with in D with pole at z, which 
solves the following Neumann boundary value problem: 

( LfiNe{-,z)) = 5z{x)I in LI, 

y^{N,{;z)} = -\dLl\-H on dLl, 

where (i^(x) is the Dirac delta function with pole at z. 
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Remark 2.24. If A € A(/i, r, k), then we have the decay estimates of the Neumann matrix as follows: 

\N^{x,z)\ < C\x - 

< \V,N,{x, z)\ + |V,iV,(x, z)\ < C\x - (2.43) 

\V,VMx,z)\<C\x-z\-^ 

for any x,z £ and x ^ z, where C depends on ix, r, k, m, d and 14. We mention that the symmetry 
condition A* = A is not necessary any longer, due to the new method developed in [2ll29j. We also refer the 
reader to [21] for the proof of (12.431) . 

3 Estimates 

Lemma 3.1 {H^ estimates). Let Ll be a Lipschitz domain. Suppose that the coefficients of Ce satisfy (11.11) 
and (II.3p . Then for any f € L^(f4;M™'^), F £ L^(14;M™') and g £ there exists a unique 

weak solution £ 14^(14;M™) to (11.51) . whenever A > K,m,d), and Aq is sufficiently large. Moreover, 

we have the uniform estimate 

11^11-5-1/2,2(9^)}) (3-1) 

where C depends only on g,m,d and 14. 

Proof. We write out the proof for the sake of completeness. First, we need to verify the boundedness and 
coercivity of B^p-, •]. The easy one is the boundedness: 

\B,[ue,4>]\<C{fi,K, m,d)\\ue\\Hi{n)\\(t)\\m{n) 

for any £ 4/^(14;M”i'). Then set (/> = Ug in (|2.5p . and a routine computation gives rise to the coercivity: 

Bs[Ue,Ue]> ^ j \VUe\‘^dx + [X- C{p.,K,m,d)] J \Ue\‘^dx > c\\Ue\\ffi^Q^^, (3.2) 

where c = min{/x/2, X — C{g, k, m, d)}. So we can choose a sufficiently large number Aq = Xo{ix, k, m, d) such 
that c = [if2 whenever A > Aq- 

The next thing is to prove (J £ = [i7^(14;M™)]*, where ^ denotes the right-hand side of 

p2.4p in the sense of 

<d,f>>=- f /" • V())"dx f F°‘(l)°‘dx+ <g,(t) > 

J J 

for any (j) £ F4^(14;M”*). It is apparent to see that 

I < 5,0 > I < c{||/|F2(o) + + II5IIb-i/2, 

and this leads to 

- C'{ll/IU2(n) + -h 115115-1/2,2(9^)}) (3-3) 

where C depends on m,d,Ll. 

Finally, due to the Lax-Milgram theorem, there exists a unique weak solution £ iL^(14;]R™') to (jl.5p . 
such that B^[us,(p] =< 5, 0 > holds for all 0 £ F4^(14;M™'). and the estimate (13.ip follows from (13.2p and 
()3.3p . The proof is done. □ 

Remark 3.2. The adjoint operator 42* has the same results as that in Lemma l3.II and the corresponding 
compatibility condition (12.6p . 
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Lemma 3.3. Let 2 < p < oo. Suppose A € VMO(M'^) satisfies (ll.ljl and (jl.2j) . Let f G 
F € L'^(Q;M"^) and g € where q = Then if F and g satisfy the compatibility 

condition f^F°‘dx+ < g°‘,l >= 0, the weak solution G to Lfiue) = div(/) + F in LI and 

duefdu^ = g — n ■ f on dLl satisfies the uniform estimate 

l|V«£||LP(n) < C{\\f\\Lp{n) + ll-^llL'?(n) + Il5lls-i/p,p(ar2)}> (3-4) 

where C depends only on g,oj,m,d,p and Ll. 

Proof. In fact, the estimate (13.4p holds for 1 < p < oo. The proof can be found in |24[ Theorem 1.1] and 
thus is not reproduced here. □ 

Remark 3.4. Before approaching the proof of Theorem ll.il we first introduce the following interpolation 
inequalities. For any u € with 2 < p < oo, and for any <5 > 0, there exists a constant Cs 

depending on 6,p,m,d and LI, such that 

ll«llLp(n) < ' 51 |Vu|| 2 ,p(o) + C 5 ||u||i 2 (j^) and ||Vu|| 2 , 9 (o) < (^||VM||ip(Q) + Cs\\L7u\\l2^^), (3.5) 

where q = Outline of proof: (i) The first inequality of (|3.5|) follows from a contradiction argument, (ii) 
The second one comes from Holder’s inequality and Young’s inequality with 6 in the case of p > while 
it is trivial for 2 < p < provided C 5 > 1 . The details are left to the reader (or see [T]). 

Proof of Theorem 11.11 The case of p = 2 has been investigated in Lemma l3.ll We first consider the 
case of 2 < p < 00 , and rewrite the (11.51) as 

( Lfiue) = div(/) + F in H 

I §z7 = 9-n- f on dLl, 

where 

/ = r + Vsu„ 

F = F - B^SIUe - (Ce + A/)Ue. 


On account of (|2.6I) . we have f^F°‘dx+ < 1 >= 0. It follows from (13.41) and (13.51) that 

l|Vn£||ip(Q) < Oj+ ||T||l9(o) + II 5 II S-i/p,p(90)} 

<C'{||/|Up(o) + l|i"llL.(n) + ll5ll B-i/p.p(an) + IkelliPfO) + 

< C'{ll/llLP(r2) + ll-^llL'JiO) + I|5 |Is-i/p,p( 90) + ll//i(0) } + <^<^11 ||iP(n), 

and this leads to 

II^£|Ivk1'P(o) < 2||V«£||x,p(q) + C'i||u£||i2(Q) 

^ C"! ll/llLP(r2) + l|-^llL-?(r2) + 115115-1/?, p(ao)}> 

where we choose S such that CS < 1/2, and note that B~^/P’P(dLl;R"^) C B~^‘^^’^(dLl;R”^) holds for 
2 < p < oo. We also use the estimate (|3.1|) in the computation above. 

In the case of 1 < p < 2, we employ the duality argument. For any h G Cq(H;M™''’*), let G 
IT^’^'(H;M™') n H^(H;M™') (p' = be the weak solution of C*(vir) = div(h) in Li and B*(v^) = 0 on 
dLi. In view of (I2.10p . we have 



VUs ■ h =< Ue,Cl{Ve) > = < Ce{Ue),Ve 


> + < > 


/ / • Vvedx + / FVedxF < g,Ve > = :< ^,Ve > ■ 
JQ Jfl 
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This coupled with 

I < i?!^e > I < ll/llLP(r2) ll^'^^ellLp'(o) + l|-^llLP(r2) ll^ellLp'(o) + C'||fl'llB-i/p.p(ao) ll^ellvyi.p'(O) 

< C{\\f\\Lp{n) + ll-^llLp(r2) + 11511 B-i/p.p(ao) } II^IIlp' (o) 


indicates 

l|Vu£||x,p(n) < C'{||/||ip(Q) + ||-F||ip(f7) + ||5llB-i/p.p(ao)}- (3.7) 

We proceed with the same idea to derive 

IkellLP(O) < C*! ||/||j;^p(Q) + ||T||x,p(f2) + ||5llB-l/p.p(ao) }; (3-8) 

and this together with (j3.7jl gives m in the case of 1 < p < 2. 

To see (j3.8p . for any H G M™'), let G M™) n {Q;W^) be the solution of Cl{w^) = H 

in fl and B*{we) = 0. Then we have UeHdx =< We > and 

I < 5^,> I < C{\\f\\LP{n) + II-^IIlp(c) + ll5lls-i/p,p(ao)}||77|| LP' (O) ■ 

The estimate (I3.8p straightforward follows, and we complete the proof. □ 

Corollary 3.5. Suppose that the coefficients of satisfy the same conditions as in Theorem W.li Let p > d 
anda = I—d/p. Assumeue G is the weak solution to (|1.5I) . where f G 

F G with q = and g G Then we have the uniform Holder estimate 

ll^e|lc0-‘^(O) < C*! ||/||ip(Q) + 11-^11^9(0) + ll5llL°°(aO)}) (3.9) 

where C depends on k, X,m,d,p and fl. 

Proof. The conclusion is immediately derived from the Sobolev embedding theorem, and we omit the proof 
here (see [Ml Corollary 3.8]). □ 

4 Lipschitz Estimates 

Lemma 4.1. Suppose that A G A{g,,T,K). Let f = (/") G F G and g G 

where p > d and a G (0,1). Then if F and g satisfy the compatibility condition f^F°^dx + 
IdQ 9 °^dS{x) = 0, the weak solution to L^{u^) = div(/) + F in Ll and due/du^ = 9 — n- f on dkl satisfies the 
uniform estimate 

I|Vu£||loo(q) < C{||/||co,a(f^) + ||T||ip(Q) + ||5llc0’'"(9O)}> (4-1) 

where C depends on fi,T, K,m,d,p,a and Ll. 

Proof. Since the estimate (14.11) has been proved in [241 Theorem 1.2] in the case of / = 0, it reduces to prove 
the estimate (14.Ih for the weak solution to the Neumman problem: 

Le(ue) = div(/) in 12, due/dvi; =—n ■ f on 912. (4.2) 

Owing to the fact that Ue in (|4.2|) can be formulated by the Neumann matrix as follows 

ufix)= f Ne{x,y)dw{f){y)dy - [ Ne{x,y)n{y) ■ f{y)dS 
Jn Jan 

= - [ ^yNe{x,y) ■ f{y)dy 
Jn 
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for any x G 11, (for convenience, we omit the upper index.) and we have 


Vue{x) = - f Vx'^yNe{x,y)[f{y) - f{x)]dy - f{x) f VyN^{x,y)dy 
J Q. J D. 

= - [ '^x'^yNe{x,y)[f{y) - fix)]dy - fk{x) f VxNe{x,y)nk{y)dS{y), 
Jn Jan 


(4.3) 


where the summation convention is used to the subscript k from 1 to d. To deal with the second term in 
the last equality of (j4.3j) . we construct the following elliptic systems 


{ Leive.k) = 0 in H, 

—-= Ukl on oil 

OUe 

for k = I,-- - jd, where Uk is the fc’th component of the unit outward normal vector to 911. In view 
of |24( Theorem 1.2], we arrive at 


< ^117111(70,7,(9^^) < C, 


(4.4) 


where we note that ||7 i||c-o,7,(9q) is bounded by a constant depending on ry, Mq. We also point out that the 
integration of Uk on 911 vanishes, which guarantees the existence of v^^k- In addition, v^^k is a matrix-valued 
function. 

Proceeding as in the proof above, the transposed matrix of k which is denoted by v* ^ can be formulated 
by 

^£,fc(a;) = / ^e{x,y)nk{y)dS{y). 

Jan 

Inserting the above expression into the formula (14.3p gives 


Vueix) = - [ Va^WyNe{x,y)[f{y) - fix)]dy - /fc(x)Vu*^.(x). 

Jn 

Then on account of (12.4311 and (|4.4I) . we derive 

l|Vu£||77c3o(f7) < C'[/]( 70 ,<T(f 7 ) -|- ||/||L°o(f7)||V7;*||7,cxj(Q) < C||/||( 70 .<T(f 7 ) (4.5) 

by noting the fact of || Vu* || 7 ,oo(q) = ||V7;£||7,oo(q), where C depends only on y,T, K,m,d,a and H. 

Combining the estimate (|4.5p and [241 Theorem 1.2] leads to the desired estimate (14.11) . The proof is 
done. □ 


Theorem 4.2. Let be given in dLll). Suppose that A G A{y,K,T), and V satisfies (II.2p and (11.41) . 
Then we obtain 

l|VTe_o||L°°(fi) < C, ||'he,o —-I||L°°(n) < C'eln(ro/e-I-2), (4.6) 

where C depends only on y,T, K,d,m and H. 


Lemma 4.3. Suppose that A G A(/i,r, k). Let he the solution to L^{ue) = 0 m H and du^/du^ = g on 
911, where 



and gij G C^(9H;M™’). Then 

C - 

|V7l£(x)l < 

for any x G H, where 6{x) = dist(x,911), and C depends only on y, K,T,d,m and H. 


(4.7) 
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Proof. The proof of this lemma can be found in |24t pp.920]. Here we can offer an easier proof based on a 
better estimate obtained in (12.431) due to the remarkable progress achieved in [2l|29]. Proceeding as in the 
proof of m Lemma 6.2], we only need to show that 

\ue{x) - Ue{y)\ <C'^ \\gij\\L°°{dn) , 

ij 


where \x — y\ < cr and r = 5{x), since the desired result straightforward comes from the interior Lipschitz 
estimates (see [3l Lemma 16]) for 


\'^Ue{x)\ < \Ue{y) - Ue{x)\'^dy) ^ . 

0\X) \J^ 

According to the solution given by the related Neumann matrices, we have 

Ue{x)-Ue{y) = / {Ns{x,z) - Ne{y,z)]g{z)dS{z) 

Jdn 

/* O 

= “N / {^i-5 - nj—){Nsix,z) - N^{y,z)}gij{z)dS{z), 

“r Jon <xzj ozi 


where we use the fact that — nj-^ is a tangential derivative in the second equality (see [211 pp.921]), 

and this leads to 


\ue{x) -Ue{y)\<2 f \V^{Neix,z) - N^{y,z)}\dS{z)'^\\gij\\L<^(^QQ) 

JdQ Y 


< 




'ijllL^idn) 




+ 


-■ ‘^'^\\9ij\\L-°{dn){h + h] 
\z-Q\<cr zf .. t J 


(4.8) 

where Q G dPl such that \x — Q\ = r, and = {z G dQ : 2^cr <\z — Q\ < 2^'''^cr}. Note that we use the 
third estimate of (I2.43P in the second inequality. 

For Ii, it is clear to see that jx — zj > r since — Q] < cr and \x — Q\ = 5{x) = r. Thus 


h = 


’\z—Q\<cr \X Z\ 


\x — y\ s , . (cr)^ ^ ^ 

jdS{z) < (cr) • — < c. 




(4.9) 


For I 2 , observing that jx — zj ~ \z — Q\ (in other words, they are comparable by triangle inequality) 
whenever z G Then we have 


00 „ 


ja: - y\ 

.v, lx — zF 
I I 


dS{z) < ^(cr) 


(2^+^cr] 


d-i 


k=0 


(2^cr) 


d — 


< c. 


(4.10) 


Plugging (14.9|) and (|4.1UI) back into (|4.8I) . the estimate (|4.7|) thus follows, and we complete the proof. □ 

Lemma 4.4. Suppose that A G A(//, k, r), and V satisfies (II.2p and (jl.4j) . Then the weak solution to 
(HID satisfies the estimate 


|VT,,o(x)| <C 


1 + 


<5(x)J 


(4.11) 


for any x G H, where C depends only on g, r, k, d, m and Q. 
Proof. Let 


Yo(^) = ^“o(^) - - ^xf{x/e). 
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By definition of and xo (see (jl.Sp and (I2.1j) ). we have = 0 in fl. If we can verify that there 

exists such that 

Kfh^ian) < Cs, (4.12) 

^ ij ^ * 


then the estimate (14.lip will follow from Lemma 14.31 directly. 

Hence the main task now is to show (|4.12l) is true. We hrst have 


dv. 


= rii 


vr-vr(«)-atf(y) 






dyj 


= nib'^oiy), 


where y = x/e. Then from Lemma 12.81 it follows that there exist such that 


We thus arrive at 




^^7 _ 
'^io — 


A 

dy.i 




and 


^“7 _ _ k;«7 
^jiO ~ ^ijO- 



(4.13) 


It remains to prove the inequality of (I4.12p . Since V satisfies (II.2p and (II.4p . it is not hard to see 
that the related corrector xo is a Holder continuous function. Therefore it follows from Lemma 12.81 that 
Ej;ieL^{Y). 

Finally, by setting g'^jix) = £E‘^^q{x/£), we complete the proof. □ 


Proof of Theorem 14.21 We proceed to prove the first estimate of (j4.6h . From the estimate (j4.1ip . it 
immediately follows that || V'I'£^o||i,°°(C 2 ) < C, whenever 5{x) > 2e. 

If 6{x) < 2e with 0 < e < 1, we will employ the blow-up method to approach this case. Let 
w‘^^{x) = H^q{£x). After a standard manipulation, we derive Li{w) = 0 in Hg and 


= ni(ex)a"/(a:)^^ = 1 =' 

^ ij ^ ^ 

where fig = {x € : ex G H}. Since clH G and g G [r, 1), the outward unit normal vector is Holder 

continuous. This together with \\'^E‘^^\\(^o,t(^y) ^ C, implies \\h\\co,r(^gfi^f^B(o, 2 )) ^ C'e. Hence from the 
Lipschitz estimate ()2.29p . it follows that (for any B{P, 1) D fie with P G dfl^, we may assume P = 0 by 
translation) 


||Vu;||Loo(B(o,i)nr2e) ^ C*! ||/i||(^o.r(gQ^piB(o,2)) + (/ \Vw\'^dx'] \ 

'' VB(0,2)nOe ^ ^ 

< Cp\£+(-f iVwl^dx) ^ 

VB(' 0 . 2 'inO. ' ’ 


where 0 < p < 1, and we use the convex property (see [THl pp.l84]) in the last inequality. Set H^fi{x) = 
(H“q(x)). By change of variable, we acquire 


I VP, 


£,0 IlL' 


dB(o,£)nn) < Cp 


|l+(/ \VH,,o\^dx) 

^ VB(0,2e)nr2 ^ 


1/PI 


(4.14) 
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It remains to estimate the second term in the right hand side above. By substituting the estimate 

Ce 


|V/f,,o(x)| < 


5{x) 


(4.15) 


into (j4.14p . we obtain 


U <CeU 

VB('0.2e')nr2 ^ VSf0.2e')nr2 


/B( 0 , 2 e)nr 2 


VS(0,2e)nr2 


<C, 


where C depends only on k, r, d, m, rj and 11. Up to now, we have accomplished the first estimate in (j4.6p . 
We now prove the remainder of (14.61) . Due to |15l Lemma 7.16], it is not hard to see 

irr . N rr . /■ , f |V77,,o(^)l 

|77,0(x) - HM\ < C{ . 


'^dz] 


(4.16) 


for any x,y G D. Let vq be the diameter of D. Clearly, we only need to estimate the first integral in the 
right hand-side, and the other follows by the same way. In view of (j4.15p . we have 


[ 

In \x-z\<i-^ JB{x,ro) \x-z\<i-^ 

dz 


< Ce 


L 


B{x,ro—£) ^ 

ro-s 


I 


^6{z) 

dSrdr 


+ C 


L 


dz 


M-i 


+ 


dB{x,r) ^(^0 ^) JrQ—e JdB(x,r) X 


B(x,ro)\B(x,ro—£) \X Z\ 

'■0 r dSrdr>, / W / ^ i 

<C'{eln(ro/e) + e}, 


(4.17) 


where we use the observation that Q C i?(x,ro) and 5{z) actually becomes tq — \x — z\, (which do not weaken 
the singularities of the integrand.) We plug (I4.17h back into (|4.16p . and have 


\He,o{x) - Hefi{y)\ <C£\n{ro/e + 2). (4.18) 

Note that defined in (II.8p is unique up to a constant. If we assume that '^£^ 0 ( 3 ^ 0 ) = / for a fixed 
point xo G D, then H^fi{xQ) = —£Xo{xo/e). So by setting y = xq in (|4.18l) . it is clear to see \H^^q{x)\ < 
Celn {ro/e + 2) for any x G D, and the desired result follows. 

As a comment, the result for e > 1 is in fact the case of the standard regularity theory (see (I2.30p ). The 
proof of the theorem is now complete. □ 

Corollary 4.5. Assume the same conditions as in Theorem 14.21 then we have 

[^£,0 - l]co,.^n) ^ C[£\n{ro/£ + 2)]^-'^ (4.19) 

for any a G (0,1), where C depends only on yL,T, K,m,d,a and D. 

Proof. In virtue of the interpolation inequality, we have 


I 1 — CT 


['^'e,0 7]^o < 4ro||'k£^o -^IIloo(q)||V'I'£_0||j;^oo(q) 


<C[eln(ro/e + 2)]i-" 

where we use (|4.6p in the last inequality. 

Corollary 4.6. Assume the same conditions as in Theorem 14.21 then we have 






< Cmaxje ’^,1}. 


Furthermore, q exists and satisfies the following estimates: 


2/3 < 


-i| 


e,0|lL°°(O) 


< 2 , 




-i| 


e,0|lL°°(O) - 


□ 


(4.20) 


whenever £ < £ 0 , where £0 = £o(f, t, n, d, m, D) is sufficiently small. 


(4.21) 
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Remark 4.7. The proof of the estimates ()4.20j) and (I4.21j) is quite similar to that given for [361 Lemma 
4.9], and is thus omitted here. 

Lemma 4.8. {A nonuniform estimates). Suppose A G A{fa,T,K), V satisfies (|1.2I) . (jl.4p . and B,c satisfy 
n . Leta€ (0,r] and p > d. Assume f G C'0’'"(O; F G LP(Li;M™) and g G M™). Then the 

weak solution to dLl]) satisfies the estimate 


llVMell^o.f < Cmaxlez 1}{ ||/||c'0.a(Q) + \\F\\lp(^q) + ||5'IIco.'^(9C)}j 
where C depends only on fi,T, k, X,p,d,m and 12. 


(4.22) 


Proof. The idea of this proof can be found in [361 Lemma 4.10], and we provide with a proof here for the 
sake of completeness. We focus on the case of 0 < e < 1, since (|4.22p follows directly from the Schauder 
estimate (|2.35p when e > 1. We hrst establish 


|VMe|| 2 ,oo(Q) < Ce”" ^{||/||co.‘^(o) + II-^IIlp(o) + 115II (an)}- 


(4.23) 


To do so, let Ve = Ue — Ue{P) — exq'^{ x/e)u'l{P) for any P G dil. By translation we may assume P = 0, 
and have the new systems: 

f C.e{ve) = div(/) + F in 12, 
yBeive) = g - n-f on 912, 

where 

/= f + £V^XO,eUe{0), 

F = F - [BeVyXO + {Ce + XI){I + £Xo{x/£))]ue{Q), 

g = g -n- VeufiO) - nA^VyXoUeiO)- 
A routine computation gives rise to 


TZ{£; 2Z2; 2A; f-,F-,g) < (77^(1; D] A] f; F] g) + C'||n£;||^oo(^) 
<C'||n,||ioo(f,)+C7{||/|bo, 

where TZ is dehned in (I2.27h . and we note that [5]c0’'^(an) ^ blco-'^ian) + C{fj,,T, K,m,d)£~^\ue{0)\. 
Then in virtue of the Lipschitz estimate (I2.32p on £ scale, we acquire 


(4.24) 


\\'^Xe\\L^(D{o,e)) A beP) + ^(^5 212; 2A; P;/; 5 ) 

< Ce®" ^[n£]co,<T(Q) + (711 tie 1 1^00 (Q) + C'{ll/llc' 0 ’‘^(n) + ll-P'llLp(r 2 ) + 11511 0 °-'^ (an)} 

ll-^llLp(n) + lbllco-'"(an)}) 

where we use (j4.24p in the second inequality, and (|3.9I) in the last one. The covering argument (the interior 
estimate is established in [361 Lemma 4.10]) together with 

Vtle = VVe + Vj/XOli£(0) 

leads to the estimate (I4.23p . 

Proceeding as in the proof above, we only consider boundary estimate, and it is not hard to derive the 
following estimate from (I2.33p : 


[V' 


^£Jco.'"(D(£)) - 


1 

121 2 


D{2e) 


+ P(e;2P;2A;P;/; 


< Ce ^[tt£]co-‘^(n) FCe ||ne||2,cxj(Q) + ||/||co.'^(n) + Ibllc0"‘^(an)} + C'll-^llLp(n) 

< ^{ll/llco-'^(n) + ll-^llLp(n) + Ibllco-'^ian)}) 
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(where we also use (j4.24jl in the second inequality and (j3.9p in the last one) and this together with the 
corresponding interior estimate (see |36l Lemma 4.10]) implies 

<C£ ^{||/||co-'^(c) + ll-^llLp(r2 ) + I|5'IIco.<^(9C)}- (4.25) 

Incorporated with (|4.23l) and (I4.25p . a simple interpolation ineqnality indicates 

11 ^ 

< 2||Vu£||2oo(Q)[VM£]^o.a(f^) < (^£2-1111/11^0,^(^2) + ||F||ip(Q) + ||fi'||co.'"(an)}) 

and then the estimate (j4.22p follows. We now complete the proof. □ 

Proof of Theorem II.21 As we mentioned before, the estimate (jl.7p immediately follows from the 
standard Lipschitz estimate (see (12.301) 1 for e > e*, and we only need to consider the following 
transformation 

uf(x) = (4.26) 

for £ < where e* = minjeo,ei,£ 2 }) and £0 is given in Corollary 14.61 and £i ,£2 will be hxed later. 


On account of (j4.26p . the Neumann problem (jl.5p can be transformed into 


L^{ve) = div(/) + F in n, 

dvp 


■ 7 ^ = 9 -n-f 


on dFl, 


where 


/“ = /“ + 

+ V^^Vv2 - BfVu2 - (c"^ + A(5"^)nf, 

=g^ -n- 

In view of (11.81) and (|2.6I) . it is not hard to verify 


[ F’^dx+ [ g’^dS{x) 

Jn Jan 


= 0 . 


Hence, applying Lemma l4.1l to (|4.27l) . we have 

l|Vu£||icx>(t 2 ) < C{\\f\\co, •'(n) + II-^IIlp(o) + ll^llcO’^'Can)}! 


(4.27) 


(4.28) 


where v will be given later. 

Apparently, our task now is to estimate the right hand side of (|4.28p . Because /, F and g much involves 
and their derivatives, based on the previous results obtained in the paper, we need to provide some 
quantity estimates related to them as follows. It follows from Theorem 11.11 and Corollary 13.51 that 


max{||M£||i4/i,p(f2), ||m£||cO.<t(q)} < (^{ II/IIL°°(f2) + II-P’IIlp(Q) + lbllL°°(ef2)}- (4.29) 

This, together with Vp = and = V('k/Q)u£ + leads to 

max{||ti£||wi.p(o), ||z^£||c'o.<^(q)} < (^{ll/lli,°°(0) + II-^IIlp(O) + llfl'llL°°(ar2)}! (4.30) 

where we use (I4.2ip in the above inequality. 

We now in the position to approach the estimates associated with F, f and g. 

In view of (j4.6p . ()4.29p and (j4.30p . we first have 


ll-^llLp(n) < (^{II/IIl°°(o) + ll-^llLp(fz) + llfl'llL°°(ao)}- 


( 4 . 31 ) 
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For the estimate of ||/||cO:‘"(n)) we rewrite as 

S/Ve = -V^efiVe + {I - 'I'£,o)Vt>e + S/u^. 

Then set u = min{r, (t}/ 2, and it follows from (14.611 . (14.191) . (I4.20p . (14.221) and (|4.30p that 

+ [I — 'I'£,o]co-‘^(c)l|V^£llL°°(n) + 11-^ “ 'I'£,o||L°°(c)[V^£]co-‘^(n) + [V^£]co.‘^(n) 

< C{£ "^ + 62 ^){ ||/||cO-'^(Q) + ll-^llLP(r2) + Il5llc0-‘^(an)} + C'l|V'f^£||L°°(n) + C£^[VVe]cO,’'(n), 
and this further indicates 

[Vu£]^o.i.(f^) <Ce’' {||/||c'0.‘^(n) + II-^IIlp(q) + llfl'llc0’‘^(9O)} + C'I|Vu£||2,oo(q), (4.32) 

whenever e < £i, where £i = min{l/(2C')^, 1} and v' = max{r, 1 — o'/2}. We thus have 


Il/llc0’‘'(f2) ^ Il/llc0’'^(n) + Ce\ii{e/rQ + 2)UVueH^oo^f^) 

+ [^£('I'£,0 - + ||K('I'£,0 — I)Ve\\cO’‘'{Q) 

<Ce^-'^'-^°\\Vvsh^^n)+C{\\f\\co, <^{£1) + ll-^llLp(n) + Il5'llc0’'^(an)}> 


(4.33) 


where we use (14.6p . (I4.19p . P4.30D and p4.32p in last inequality, and choose a small positive number lq such 
that e''° ln(ro/e + 2) < 1. 

We now proceed to show the quantity of ||3||c0"‘^(an)) it is not hard to see that 


ll^llco.-^Can) < llfi'llco.'^(9n) + C'll'f^£llcO''^(c) 

< C{\\f\\c°’<^{n) + II-^IIlp(c) + Ibllc0’'^(ar2)}- 

Collecting (I4.28p . (I4.3ip . (I4.33P and (|4.34l) yields 

l|V^£lli,°°(r2) < Ce^ + Cj ||/||c-o.<T(Q) + + ||fl'llco-‘^(9r2)}! 

and this implies 

l|VU£||ioo(Q) < C'{||/||co, '^(n) + II-^IIlp(o) + Il5llc0’'^(an)}) 


(4.34) 


(4.35) 


whenever e < £ 2 ) where £2 = min{l/(2C') ^ 1}. Recalling the transformation p4.26p . it is straightfor¬ 
ward to derive the estimate ()1.7p from (|4.35l) . Up to now we have proved this theorem. □ 


5 Convergence rates 

Lemma 5.1. Suppose that Ue,uo € F1^(U;M"*) satisfy Cei^e) = ^o{uo) in fl, and Bs{us) = Bo{uo) on 50. 
Let 


= u^-u^-e'^ xi'^{x/e)ipl, 
where ip = {pj.) G Then we have 


k=0 


d 


= -jF R? + K" - <€] 


dx 


du\ 


0 h 
--pf 


dx, 


+ [U'’ - v°f] K - *> 0 ] -4rr + 3t) 


rdU, 


+ - bZ] [g - ^f] + - cf ] H - Z] -e{M-+ 


(5.1) 
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and 


[8.(».)]” = n,{lCt + + [!/“'’ - V.f ] [,.$ - f.?] - elf}. (5.2) 


where Ui denotes the i’th component of the unit normal vector to dQ, and 

sn 

&U: 


+ Jt = T.-7t^l K? = ^iq:p], 


k=0 


k=0 


k=0 


j=0 


d 


•^“ = E[ 

k=0 


ddV _ R fl .,1 d 


(5.3) 


dy. ^,e Q^. 

Proof. In view of Cs{u^) = Cq{uq), we have 


+ f-Pl}, V” = [cf + A5“'>] 5: xfjri. !/ = x/e. 


A:=0 


A:=0 


where 


= Eo(no) - ^^(uo) - Ee (e ^ Xk,ePk) 

idu^ 


(5.4) 


[£o(xo)]" = -£{ihf^ + K '‘<}++ r" + 


dxi 


and 




(5.5) 


- .B. 


- ^Icf + X4“'’lxf>2 


(5.6) 


with summation convention applied to k from 0 to d. By substitnting ()5.5I) and ()5.6p into (15.41) . we obtain 
\ r f Al “ ^ f ''OiP^UQ ajj 9^0 0/3 ^Xi 7 , Tyo/3 /3 txo/ 3 /3 o/3 ^Xo 7 

[^'('"'>1 = ^ < - ''‘f< - 

'1 + po/j ^^0 _ _ paS ^xh 1 , _ poA ^Xq^ 


-eZf 


7 

V^o 


(5.7) 


- =<Exr 7 {v 2 } - E|ef + X 5 “'’] Exf>/, 

k=0 * fc =0 

where / = 1,..., d. By definition of (12.121) and (I2.13p . we have 

^a/S^^O ap ^^0 0/3 ^X/ 7 ,07 7 , r^o/3 o/3 1 P 

“« ^ = "«>? + [“« - “«.J l^-Ri 

vrV - K°V - = tpiPi + |l>“'’ - VPV„ - P>„] 


(5.8) 


'i “0 i,£ “0 ~ Qy . 


and 


i®“o „oeS"o pae^xfx p 




2,£ 

/7 


dxi 




- cf 4 - + [c-f - cfU - 


(5.9) 
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Then collecting (I5.7p . ()5.8p and (I5.9p gives 


[£.(».)]“ = -^{ E C*’2 + lay - <1 [|^ - *^] + IT" - v',?ll»g - - £i."} 




fc =0 


0/3 730 / 3 ] 


2 ,£ • 


+ [^7 " - B\, 


du^ 


d 


Vdxi 

- e[cf + XS‘-f‘]j2xll‘Pl 




k=0 


dxi 


A:=0 


A:=0 


This together with 


Ewk*’ 2 -E|:{^ 4 -E 


k=0 


k=0 


dxi I dy, 


k=0 


d 

dyi dx, 


{^1} 


leads to the expression (15.11) . 

We now turn to prove (15.2p . Note that Be{us) = 13o{uo) on dQ, and we have 


Be{We) = Bo{uo) - Be{uo) - Be{e'^Xk,e^k) 


where 


and 


[Bo{uo )]" = [Beim )]" = nW“/4 + 


k=0 


^ . T /“/ 3,,/3 1 ^ *^^0 


dx 


ri 

[Be{exk,s^k)r = + 

for A: = 0,1,... , d. Then inserting the expression (15.111) and (15.121) into (j5.10p indicates 


d 


(5.10) 

(5.11) 

(5.12) 


Th / \1 “ /t/o/3 /3 t/o/3 /3 0/3 *9X0 7 , ~o/3*9^o ofi ^^0 o/3 “^X; 7 .y, 

[8,(»,)] = „„ - 1/.., „„ - + »y ^ - »y,,^ - 


?}, 


where summation convention is used to I from 1 to d. This together with (|5.8p finally gives the expression 
(15.2p . and we complete the proof. □ 

Lemma 5.2. Let and uq be given in Lemma [5Tl and is defined in (|1.12l) and satisfies (|5.1I) and (|5.2p . 
Then we have 

ll^£‘?llL2(0\E2e) + 11^3^0 — fi\\L^{n) + II^^O “ To\\L^{n) 

^ (5.13) 


+ £\\h^V(j)\\i2f^Q^ +e||/i£ 


+ Ce\\he4>\\L'2[dn), 


where (p = ((/?i,--- ,pd), 4> = h represents the periodic function depending on some of the 

periodic functions shown before such as the coefficients of the correctors {xk}k=Q, and auxiliary functions 
{bik ) Ejik ) ^k}k=0 • 


Proof. Let We = where Wep and satisfy 


Ceiw£,i) = Bfiwe) + div(/C) in 12, 
Be{We,l) = Be{We)-n-{lC) on (912, 


(5.14) 

























and 
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C-e{We, 2 ) = -div(/C) 

Be{We, 2 ) =n- 1C 

respectively, where 1C = (/Cf) is defined in (j5.3p . 

For the hrst equation in (j5.14p . it follows from (j.S.ip that 

— C'{||Vno — + ||no — <fo\\L‘^{yi) + A\he'^$\\L‘^{yi) + A\he$\\L‘^{yi) + A\he$\\L‘^[dQ.)] ^ (5.16) 

where h depends on the coefficients of the correctors {xk}k=Q and auxiliary functions {'C^'&k}k=o- 
We now focus on (I5.15p . In view of ()2.4I) . we have 


in n, 
on dCl, 


(5.15) 


^e{we, 2 iv\ = / 1C ■ Vvdx =: R{v). 

Jn 


According to Lemma [2181 R{v) in (|5.17l) satisfies 

= " X ^ ^ - £ X E 

k —0 k —0 

=: Ri{v) - R 2 {v). 

Note that due to the antisymmetry of Ejik with respect to i,j, we obtain 


= e 


= e 




d'^i 


dxidxj 


dx 


where ?/)£ € (7^(11) satisfies (I2.36p . and A; = 0,1,... , d. Moreover, we have 

= ■ X + X<' ■ ^ 

and this indicates 


dtfl 
^ dxj dx. 


(5.17) 


-dx, 


|d2l(n)| < {||L^e(X||L2(f2\22e) + ll^e'?llL2(0\S2e) +^l|-^£'^<?llL2(n\S2e)}ll'^^llL2(0\S2e)' (5.18) 

Meanwhile we arrive at 

\R2{v)\ < s\\Ee\/$\\L2(^a)\\Vv\\L2(^^y (5.19) 

Let V = We ,2 in the expression (I5.17p . Combining (15.171) . (I5.18p . (I5.19P and (|3.2p . we acquire 

11 ^£,2 11//I (fi) < C'{l|-E^£<?llL2(n\S2e) + ll^£<^llL2(0\S2e) + (U)} ■ (5.20) 

Then the conclusion of this lemma immediately comes from (|5.16l) and (I5.20p , and we complete the proof. □ 
As a reminder, Se^Se^tpr is defined in (I2.2np . (j2.24p and (j2.36p . respectively. 
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Lemma 5.3. Let Q be a bounded domain. Suppose that the coefficients of satisfy (|l.lh . (jl.2p and 
(USD. Assume that u^, uq are the weak solutions to TO with F e and g G Then 

we have 

llue -Uo- SXO^eSeilpAeUo) “ £:Xk,6Se{'lp4e'^ kUo)\\< <^£^11^11^2(0) + |bllsi/2,2(an) } : (5-21) 

where C depends only on g,K,m,d and Ll. 

Proof. By choosing ipQ = Ssif’AeUo) and tpk = Ssi'ipie'^kUo) in (I1.12p . we let 

We = Ue - Uo - eXO,eSe{'4’AeUo) “ eXk^eSeif’Ae^kUo)■ 

Then it follows from Lemma 15.21 that 

I|w^£||h 1(!^) ^ ll'^O “ Se{i^AeUo)\\L-^(n) + l|Vno - Se{'tpAe'^Uo)\\L2(n) + £\\heSe{'tpA6Uo)\\L2(n) 

^ (5.22) 

+ e\\heSe{lpA 6 '^Uo)\\L 2 {n) + e\\heVSe{lpA 6 Uo)\\L 2 (n) +e||^£V5£(V'4£Vuo)||L2(Q) L 

We note that SeitpAeUo) and SeitpAe^uo) is supported in Ss^. To complete the proof, we need the following 
estimates. 

Due to (I2.22jl . we have 

11^0 - Se{'lpAeUo)\\L^{n) < 11(1 “ V’4e)';^0 ||l2(q) + Wf^AeUa “ <S'e(V'4eWo) ||l2(q) 

< lko||L2(f^\S8e) + C'£||V('04e'no)||2,2(f7) (5.23) 

< C'lko||L2(n\S8e) + C'^ll^^o||L2(S4g) 

and 

||Vtto - Se{ipAs'^uo)\\L2(n) < C'I|Vmo||l2(q\s8j) + C'£||V^uo||i 2 (S 4 ^). (5.24) 

From Lemma l2.ini it follows that 


\\heSe{'lpAeUo)\\L^n) + Whe'^Se{'lpAeUo)\\L^{n) < C\\lpieUo\\L^{n) + C\\V{lpieUo)\\L^n) 

< C'lko||L 2 (S 4 j) + Ce ^||Mo||L2(Q\28e) + C'II^^o||l 2 (S 4 ^) 

and 


||^£'S'e(V’4eVMo)||L2(Q) + \\heVSe{'lpAe'^Uo)\\L 2 (n) < C'llV mo||l 2 (s 4 ^) + Ce ||VMo||L2(Q\s8e) + C'||V tto||L2(S4e)- 


Plugging the estimates (15.231) . (15.241) . (15.251) and (15.261) into (I5.22p . we obtain 

ll'*ne||_H-i(Q) < C'lko||Hi(0\Ege) + Ce\\uo\\H2{T,4e)- 


(5.26) 


(5.27) 


Set g = {gi, ■ ■ ■ , gu) G Cq (M*^; M'^) be a vector field such that {g,n) > c > 0 on dQ, where n denotes the 
outward unit normal vector to dLl. Since the divergence theorem, we have 

c[ {\^uo\^ + \uo\^)dS< [ <g,n>{\Vuo\^ + \uo\^)dS 

Jan Jan 

= / div(^j)(|Vuo|^ + |uo|^)dx + 2 / gly^uo^uo +'VuoUo)dx (5.28) 

Ju Jfa 

< C{\\uq\\\2(q-^ + ||Vuo||i2(Q) + ||V^Uo||L2(f2)} < C'||'Wo|lH2(f2) 

By the same argument, it is not hard to see 

( [ (iVnoP + k|2)dSt) < C||no||^2(s,) < 

'J St 
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holds for any t £ (0,8e), where St and Sf are defined in Remark 12.191 and C does not depend on t. Then 
by (12.401) . we obtain 


l|wo||/fl(0\S8e) 


( / (|VnoP + |wP)dx') 

^ Jn\j:8e ' 

(/ {\Vuof+ \uf)dStdt^ ^ < Ce^ 11^011^-2(5^). 


(5.29) 


This together with (I5.27|) gives 

< <^£^1^011^2(^7) < <^£2 {||-^IIl2(o) + ||5'll_Bi/2,2(af7)}. 

We note that the assumption of G is only used in the last inequality, and we are done. □ 

Remark 5.4. We replace in the proof of Lemma 15.31 into 

We = Ue - Uq - eXO,eSe{uo) - £Xk,eSe{^kUo), (5.30) 

where ho is the extension of uq such that ho = uq on and ho G M™) satisfying ||ho||// 2 (]Rd) < 

C'||no||H2(o)- Also, Se is given in the sense of (I2.24p . It follows from (I5.13P that 

ll^e||/7i(C) < C*! ll^£‘S'e(ho)||L2(r2\S2e) + ll^£'5£(^'“o) ||L2(f2\s2e) + 11^0 “ '5e(ho) ||L2(f^) 

+ ||Vho — 5£(Vho)||L2(Q) + e||/le5e(V^ho)||L2(Q) + e||/l£5e(Vho)||L2(Q) 

+ £||^e'S'£;(ho)||i2(f7) + e||/i£5£(ho)||i,2(gn) + £||^£‘5£(Vho)||L2(9o)|, 

where we need the counterparts of (12.211) and (I2.22p for .^^(or see Remark 12.121) . The different task from 
the proof of Lemma 15^ is to estimate the terms of heSe{u) and heSe{Vu) in L^(12 \ T, 2 e) and L^{dQ) by 
Lemma 12.141 Then we can also arrive at 

{\\^\\L^{n) + 11511^1/2,2(90)} (5.31) 

without any difficulty, which actually gives the counterpart of [331 Theorem 4.2] in our case. 

Lemma 5.5. Let Q be a bounded domain. Suppose that the coefficients of Ce satisfy (|l.ip . (|1.2p and 
do]). Assume that Ue, uq are the weak solutions to dO]) with F G L2(12;M”^) and g G Then 

we have 


ll^e - Uq- eXo,£*S's(V’4eho) - £Xfc,£'S'£(7/^4£ Vfcho) || ^2(o) < <^£1 I|-F’||l 2(0) + ||ff 11^1/2,2(90) }, (5-32) 

where uq G L1^(M‘^;M™') is the extension ofuQ, and C depends only on iJ.,K,m,d and Q. 

Proof. We will employ duality argument. For any <1> G L^(11;M™’), let (pe and pQ satisfy 

I ^ in 11, f Lq((Po) = 4> in H, 

[B*ffife) = 0 on 911, }0*(,)6o) = O on 911, ’ 

respectively. Due to ()5.2ip . we can have 

\\4>e - 4>o - £Xo{x/£)Se{'ifioe4>o) “ £Xfc(a:/£)S'£('0iO£Vfc(()o)||j;^i( 0 ) < ^£^11^111,2(0), (5.34) 


where xl nre correctors associated with £*, and satisfy the same estimate as Xfc do, where k = 0, - ■ ■ ,d. 
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Let 


We = Ue - Uq - £XO,£<S'e(V'4eWo) “ £Xk,eSe{i’Ae'^kUo), (5.35) 

where uq G is the extended function of uq such that ho = uq in and ||ho||j:f 2 (]Rd) < C'||tio||// 2 (f^). 

According to the formula (|2.11l) and Lemma l5.11 we obtain 


[ f ■ '^4>edx+ [ 

J Q J Q 


/ We^dx = I f-S/4>edx+ I F(j)edx, 

J Q, J Q 


(5.36) 


where 


f = IC-e{I + J) + [A- A,] [Vuo - 5 e(V’%Vho)] + [F - F,] [uo - Sei^l^4em)] , 

F = -e{M+Af) + [B - Be] [Vuq - Sei'lpie'^Uo)] + [c- Ce] [uo - Se{lp4eUo)] ■ 

Below we do some calculations in more details. By the fact that uq is the extension of uq, we have 

Vuo - Sei'llj4e'^Uo) = Vuq - ^^(Vho) + 5'£((1 - V'4e)Vho) “ (1 “ V'4e)V'Uo + (1 - 1p4e)'^U0 in 12, 
and then 


<C<^ y^|(l-V’4e)Vuo||V(/.,|dx 


[ [A- Ae][Vno - Seillj4e'^UQ)] ■ V (pedx 

Jn 

+ [ |Vho-5e(Vho)||V0,|dx+ [ |(l-^4s)Vho-5,((l-V'4.)Vho)||V0,|dxi 

Jn Jn J 

< C’|||Vuo||L2(Qys8e)l|V</>e||L2(n\S8e) + ll^ho - *S'e(Vho) 11^2(Rd) ||V(/>£ 11^2(5.37) 

+ 11(1 - V'4e)Vho - 54(1 - V’4e)Vho)||i,2(Rd)||V(/i£||i2(Q\s9e)} 

< C'||Vuo||L2(n\S9e)l|V^e||L2(Qys9e) + Ce\\V'^Uo\\L2{M.d)\\V(j)e\\L2(n) 

+ C'£||V[(1 “ V'4e)Vho]||i2(Rd)||V(()£||i2(Qys9e) 

< C'||Vuo||L2(Q\Sg^)||V(/)e||j;^2(Qysg^) + CeUuo ||h 2(O) II V'/’e lll,2(n), 

where we use Cauchy inequality in the second inequality and the observation that 5e((l —V' 4 e)Vho) restricted 
to n is supported in 12 \ Bg^. In the third one, we employ the estimate (I2.22p . and the last one follows from 


||V[(1 - V’4£)Vho]||i;,2(Rd) < Ce ||Vho||L 2 (o\E 9 e) + dliv ho||i;, 2 (Rd) < Ce ||Vno||L 2 (o\E 9 e) + C'lko||H 2 (C!)- 

Similarly, by noting uq - Sei'ilJ4eUo) = uq - Seim) + 5e((l - ip4e)uo) - (1 - '04£)ho + (1 - '04e)Mo in 12, we 
obtain 

/ [V - Ve][uo - Sei'llj4EUo)] ■ V (j)edx < C||no ||L2(f2\s9e) II We llL2(Q\S9e) + C'^II^O Ib2(f^) || V(/>£ ||i2(f2). (5.38) 
Jn 

As demonstrated before, it is not hard to derive 

[B - Be] [Vuq - Sei'lp4e'^m)]4'edx < C||Vuo||i 2 (Qys 9 e) ll0elll,2(n\S9e) + C'sll^O 11^2(0) ||<(>£ ||L2(n), 

J (5.39) 

/ [c-Ce] [uo - Sei^p4eUo)](l)edx < C||uo ||L2(f2ys9e) ll^e llL2(n\E9e) + C'^II^O 11^2(0) I|0£ 1^2(0) ■ 

In 

Now, for simplicity of presentation, let h'^jk,e^x) = with /c = 0, • • • , d. Set V* = ^, , 

and the notation of “VoUq” means uq itself. Note that 

/ K^k,e'^,{Seiid4eykul)]Viredx= [ h“4,54V’4e V^h^) + [ h,7,^,54V,u/>4e Vfch^l) 

(5.40) 
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and then we have 

< C'(||V^uo||i 2 (Rd) + ||Vuo||i 2 (Rd)) ||V(/)£||i 2 (Q) + Ce ^||tto|| Hl(n\S9e)l|V</>£||L2(QySg^) 

< C\Wo\\H^{n)\\'^4’£\\L^{n) + C'e~^lko||Hi(o\E9e)l|V</>£||L2(Qysg^), 

(5.41) 

In the first inequality, we mention that S^(VjipieS/kU^) restricted to 11 is actually supported in H \ Sg^. In 
the second one, we use Lemma l2.ini In the last one, we note that uq is the extension of ug- 

If we set with k = 0, • • • ,d, then the following estimate immediately holds by Lemma 

Eini 


Jn 


— ll^rfe^£‘S'e(V'4£Vfc?2(()||i2(Rd)||Vi(/)"||i2(Q) < C'||uo||Hi(o)||V(/)£||i2(Q). 


This together with the estimate (|5.4ip gives 


I ■ \/(j)i;dx 


< C\\uo\\H‘2{n)\\'^4>£\\L'2{n) + Ce ||^io||Hi(n\E9e)ll^'/’£llL2(n\s9e) 


(5.42) 


by recalling the expression of I in (j5.3[) . By the same token, it is not hard to acquire 


J ■ Wcpedx 


M-Cpedx 


< C\\uo\\m{n)\\'^(pe\\L‘2{n), 

< C'lko||_ff 2 (f 2 )||<^e||L 2 (r 2 ) +Ce ^l|Mo||_ffl(f 2 \S 9 e)ll</’£||L 2 (f 2 \S 9 e)) 

/ Afcj)^dx <C\\uo\\H^n)Ue\\L^in)- 
Jn 

The rest thing is to study the term of JqIC ■ Vc^^dx. In view of ()5.17p . we have 


(5.43) 


/ /C • VcP.dx = R{(t>e) = Rl{(ke) - Rim- 
Jn 

It follows from the estimate ^5.181) that Ri{(j)£) = 0, because Se^ipAeUo) and Se^ipAe^uo) are snpported in Egg. 
The calculations of estimating i? 2 (</*£) are quite similar to (I5.40D and (15.411) . and therefore some explanations 
are omitted. 

R^m = e J^Efl^\/j{SmAe^km}^i€dx 

< C'll^^o||//l(n\S9e)l|V<('£|lL2(Q\S9e) + C e\\uo\\jj2 (^q-^WV 4>s\\L'^ (n) ■ 


This implies 


JC ■ Vcjiedx < C'||tto||Hi(!^\S9e)ll^'/’£llL2(c\S9e) + Ce\\uo\\fj2(ji^\\V(t)e\\L'2{n)- (5.44) 

Inserting the estimates p5.37p , (I5.38P , p5.39p , p5.42p , ^5.431) and ^5.441) into (I5.36P , we derive 
J^We^dx < C'|||no||Hl(o\S9e)ll</’£llHl(*^\S9£) + ^lko||//2(Q)||(/)£||j^l(f2)|, 


(5.45) 
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To improve the order of the convergence rate, the next task is to replace </>£ into the first order corrector 
in the first term of the right-hand side of (15.451) . where — 0o — eXo,£‘5'e(^/’iO£</>o) e*S'£(V’iO£Vj(/>o). 

Observing that and 5e(?/)ioeVj(/>o) are supported in Sg^, we arrive at 


ll</’e||Hi(!^\S9£) - ll^£lli/i(*^\S9e) + ll</’o||_ffi(r2\S9e) 

Thus we plug (j5.46p back into (15.451) and obtain 


(5.46) 


/ 

Jn 


We^dx 


< 


C|lko||Hl(f^\S9£)ll'i^olli^Tf2\S9e) + ll^0||Hl(f1\S9e)ll?£lli?l(f2) + ^ II ^0 11 2 (f2) ||'^£ || (f^) } 


< c 


:)ll'/’0||/fl(f2\S9e) + 1^011^1(^^X296)11^^1^2(0) + £ 


(5.47) 


where we use the estimate (j5.34j) and Lemma 13.11 in the last inequality. 

The problem reduces to estimate the layer quantity. On account of the estimate (15.291) . in fact we have 


lko||//i(o\S96) - lko||H2(o), ll</’o||i/i(o\S96) - WMn^in)- 


Thus it follows that 


/ 

Jn 


Wp^dx 




and we consequently obtain the desired estimate (I5.32p by using duality and estimates (where the 
assumption of dQ G has been used). The proof is complete. □ 

Lemma 5.6. Let Q be a hounded Lipschitz domain. Suppose that the coefficients of Cp satisfy (jl.ll) . (11.21) 
and (ILSI), and A additionally satisfies 64 = ^4*. Let Up,UQ be the weak solutions to (II.9p with F G LP(f4; M™) 
and g G L‘^{dQ]W^), where p = Then 

\\ue - no - exo,eSl{fi 2 eUQ) - eXfc,£'S'£ ('02£VfcUo)||< C'e^||T||ip(s^) -h ||fif||L 2 (ao)}, (5.48) 


where C depends only on g,K,m,d and fl. 


Proof. We note that UtOellHTo) is exactly the left-hand side of (I5.48P by setting = S‘l{'fi 2 e'^o) and = 
Sp{ip 2 e^kUo) in (11.121) . Then it follows from (15.131) that 


llw^elli/TO) < C'{l|Vuo - S^{Tp26'^Uo)\\L^n) + Iko - S‘^{fi2eUo)\\mn) + 4K'^S‘^i'^2eUo)\\L^n) 
-h £\\heVS'^{fi2e'^Uo)\\mn) + e||^£5'£ (V'2£Mo) ||l2(0) + e\\hpS‘^{fi2e'^Uo)\\L2^n)}- 
Before proceeding further, let us do some calculations: 

||Vuo — Sp{'lf 2 £'^Uo)\\L^(n) < 11(1 — V’2£)Vuo||L2(f7) -I- ||V’ 2 £Vuo — 5'£('02£Vuo)llL2(f7) 

+ IIS's [V'2 £Vuo - 5s(t/'2£Vuo)] ||l2(q) 

< l|V^o||L2(n\S46) + C'||'i/’2£Vt(o - S£('(/’2 £Vuo)||l2(q), 

\\heVS^p{fi2eVuo)h2in) < C|| V5s(V'2£Vno) ||i2(f,) 

< Cje ^||VUo||l,2(n\S4e) + l|S£(V’2£V^t(o)||L2(f2)}, 

and 

l|/i£S£('02£Vuo)||L2(Q) < C'||5 £(t/>2sVuo)||l2(0) 


(5.49) 


(5.50) 


(5.51) 

(5.52) 
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where we mainly use the estimate (12.21jl in the second inequality of (j5.50j) . and in the first inequality of 
(15.511) . as well as in (15.521) . After a similar computation, we have 


(5.53) 


\\heS‘^{'lp2£Uo)\\L2(^Q) < C\\Se{'llj2eUo)\\L^{n), 

11^0 - S^{'ilj2eUo)\\L^n) < lko||L2(0\S4e) + C\\'lp2eUo - Ssi'tp2eUo)\\L2(n), 
\\he'\/S‘^{lp2eUo)\\L^{n) < <^1“ lko||L2(0\S4e) + ll'S'£('02eV'Uo) |. 

By substituting ()5.50p . (|5.5ip . (15.521) and (15.531) into (I5.49p . we find 

- ^{ll^ollfl'i(0\E4e) + IIV’2eVuo - 5'£('02eVno) || 1,2(0) + ||V’2£^^0 “ 'S'e('i/’2e^to) IIl2(0) | 

+ C'e|||S'£(V'2eV^Uo)||L2(o) + \\Se{'4>2e'^Uo)\\L2{n) + l|5'£(V'2e'Wo) ||l,2(o) |. 

We now handle ||^io||J^l(o\S 4 e) the right hand side of (I5.54p . First of all, we rewrite (Hq) of (|1.9p as 
Lq{uo) = F + {V — uq — (c + XI)uq in 12, du^/duQ = g — n-V uq on 912, (5.55) 

where Lq = —div(AV) and djdv^ = n ■ AV. Let uq = v + p, and v, p satisfy 

in 12, 


(5.54) 


(Di) Lq{v) = F in 


(D 2 ) 


Lq{p) = 0 
dp 


dv 


= g-n-Vuo-^ on 912, 
duo dvQ 


(5.56) 


respectively. Note that F = F F {V — B)Vuo — (c + \I)uq in 12, and F = 0 in \ 12. 

For (Di), let Fq denote the fundamental solution of Lq, and then we have v = Tq * F in Moreover, 

we have the estimate 

\Hy)\ 




\x - y\<^ 

and by Hardy-Littlewood-Sobolev inequality on fractional integration (see |181 Theorem 7.25]), 

I|V»||,.„., < CIIFII „ < C{||F|| „ + llab^m). P-v) 


'■ L^(n) ^ '■ L^(r2) 

Also, it follows from the Calderon-Zygmund theorem (see [181 Theorem 7.22]) that 

llV^ull 2d <C|]F|] 2d <C{||F|| 2d + Ikolbifml < CjllT^II 2d + ||5 ||l2(9o)}. (5.58) 

L^(Rd) “ L3+T(Rd) “ "lH+T(0) “ L" "lH+T(o) " 


Note that we use (|3.ip in the last inequalities of (I5.57P and (|5.58p with the fact of L^(911; M”*) C B ^/^’^(912; M”^). 
Then due to Sobolev inequality, we have 


^ <^{11^11, + lbllL2(90)}, 


L3+2(0) 


||Vu|| 2d <C|]V2u|| 2d <C{|]F|] 2d + |bllL2(ao)}- 

1,3=0 (Rd) - "l,3+T(Rd) - Lll ^3+1(0) ^ 


(5.59) 


Hence, by suitable modification to the proof of the estimate p5.28p . we acquire 

c[ (|Vup + |up)d5< / < Q,n> {\Vv\‘^+ \v\‘^)dS 

Jan Jan 

= / div(£>)(]Vu]^ + ]u]^)(ix + 2 / giv'^vVv+ Vvv)dx 

Jq J o 

<c{\\vf +||Vu||2,,^.+ ||Vu|| ^ llV^ull ^ } 

T.X=^(0\ Hi (ip 


(5.60) 


L3=0(0) 
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where we use Holder’s inequality and Young’s inequality in the second inequality, and the estimates (j5.57ji . 
(15.581) and (I5.59P in the last one. As explained for the estimate (15.291) . it is not hard to see 

{i f (|Vu|2 + |u|2)dx|" <C7{||F|| ^ +||5||i2(a^)}. (5.61) 

Jn\T.4^ J Ld+l(n) 

We now turn to study (D 2 ). From the Rellich identity, the estimates ()3.ip and (I5.60p (see Remark 15.7p . 
it follows that 

ll(Vp)*||L2(gn) < C{\\g\\L2(^QQ-j + (5.62) 

According to Lemma 12.231 we have 

\\^{p)\\L^{dn) < C'lbll/fi(n\Ec(,) ^ C\\p\\Hi(n) < C*!lbllL2(an) + 

where we use (j3.1l) in the last inequality, and C depends on d, m, cq and the character of H. This together 
with p5.62p leads to 

{7 [ (|VpP + \p\^)dxY < C{||(Vp)*||i.(ao) + \\M{p)h2^9n)} 

(5.63) 

<C{||9|b,8n, + ||F7^^^7 

where we use the estimate (12.391) (for p = 2 and r = 4e) in the first inequality. 

Hence combining (15.611) and (I5.63p . we have 

{[ (|VuoP + |uoP)dx|" < C'e^||ff||i2(af2) + lli^ll M }. (5.64) 

We now estimate \\^p 2 e^uo — Ss{ip 2 e^uo)\\i 2 ^Q^ and \\'ip 2 eUo — Ss{ip 2 eUo)\\L^^Q) in the right hand side 
of (j5.54p . Since V’ 2 £Vuo ~ Ss{ip 2 e^uo) is supported in it is equivalent to estimating \\'ip 2 e^uo — 
>S'e(V’2£VMo)||L2(Ee)- It immediately follows from (j2.22p . (j2.23p and (I5.61j) that 


||V’2£Vuo - S'£(V’2£Vuo)||l2(s7 

< ||Vu - 5 ,(Vu)||^2(s,) + ||(V’2£ - l)Vu||i2(s7 + \\Se{{^2e " l)Vu) ||^2(s7 + ||V’2£Vp 

< 1|VU - 5£(Vu)||^2(Kd) + \\{'4>2e - l)Vu||i 2 (f^) + ||5'e;((V’2£ “ l)Vu) ||i 2 (s 7 + ||V’ 2 £Vp- 

< C£2 + C'||Vu||j;^2(Q\S4e) + C£\\V{Vpil;2e)\\L^{n) 

< C£^{\\F\\^^^^^ + ||5||L2(gn)} + <^{£11^^^11^2(225) + II^Plli2(n\S4e)}- 

Note that p satisfies (D 2 ), from the interior estimate for Lq (see Remark 15.8p . we have 


-5£(V’2£Vp)||^2(25) 

5’£(V’2£Vp)||i2(f2) 


(5.65) 


vV(x)| < 


C 

5{x) 


B{x, 5 {x)/ 8 ) 


Npiy)\"dy) 


1 

2 


This gives 


'S2e 


<C [ / 


S2e JB{x, 5 {x)/ 8 ) 


|Vp(l/)P 

[5(X)]2 


■dydx < C 


n 


|(Vp)*(xO|^ 


f2 


■dSfdt. 


(5.66) 


Note that 6(x) « t, and x' G <S' 5 (x )/4 such that \Vp{y)\ < {Vp)*{x') for any y € B(x, d(x)/8). By using the 
observation that ||(V/?)*||j;^ 2 ( 5 j) < ||(V/9)*||£2(gQ) holds for all t G [0,oo), it follows from (I5.62p and (I5.66P 
that 


^^PllL2(E2e) - PTWh’^ido.) < Ce ^{\\g\\L'2{dn) + 


(5.67) 
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Now, collecting (I5.63j) . (j5.65jl and (j5.67p leads to 

||V'2£Vmo - Se{'lp2e'^Uo)\\L2(^Q) < Ce^\\g\\L2^g^) + (5.68) 

By the same procedure as above, we have 


\\i’2eUo - Se{'llj26Uo)\\L^(n) < <^£11 V('02e'Uo) ||l 2(0) < C'lko llL2(n\E4e) + C £\\\/UqW ^2 


(5.69) 


where we use (j2.22p in the first inequality, and (13.ip and p5.64p in the last one. 

To accomplish the proof, we still need the following estimates: 

||5,(^2.V2uo)||i2(f,) < ||5,(V'2sV\)||i2(f2) + ||S,(V'2 sVV)IIl2(0) 

<Ce-k^P2sS/M^^ +C\\V^p\\lH^,.) 


(5.70) 


<C£ 2{||5||^2(aQ) + 

where we use p2.23p in the second inequality, and p5.58p . p5.67p in the last one. Also, 
\\Sei'4’2e^Uo)\\L^{n) + ||5'£('02e'«o) ||L2(n) < C'll«o||//i(n) < <^{1151^2(9^) + 


(5.71) 


Finally, combining (15.541) . (15.641) . (I5.68F (15.691) . (15.701) and (I5.7ip . we obtain the estimate (I5.48P and the 
proof is done. □ 


Remark 5.7. To see the estimate (|5.62l) . we hrst show the Rellich identity for the equation: Lq{u) = F m 
n and dujdvQ = g on dO., where F € L^(n;M™') and g G L^(cin;M™') with F°‘dx + fg^ g°‘dS = 0. 


/ 

Jan 


^ f pQk du°‘ ^ ^ 

Jn 


dx. dXi 


I 

Jn 


(9u" 

Q dxj dxi ' “io dxk 

-of 

Jq dxk dxj dx 


F'^Qkdx 


-dx 


f dvP^ 

2 / rj^QkdS, 
Jan 


(5.72) 


where ^ is a C'q(M'^;M'^) vector held similarly dehned as in (I5.28p . and we use the assumption of A = A*. 
Coupled with the estimate (see (12.311) for p = 2 and / = 0), it is not hard to see that 


l|Vu||i2(9Q) < C{\\F\\j^2(^Q) + ||3llL2(aQ)}. 


(5.73) 


Then, we apply the estimate (|5.73l) to the solution p to (D 2 ) in (|5.56p (where F = 0 and g = g — ti-Vuq — 
dv/dv), and acquire 

II^PllLRan) < C'{lbllL2(ao) + ||'Wo||L2(afi) + l|V^^llL2(af2)} 

< C{\\g\\L'2(an) + IkoliHRCJ) + l|Vu||2,2(gf2)} < <^{1151^2(9^2) + 

where we use the trace theorem (see m Lemma 5.17]) in second inequality, and the estimate (|3.ip for uq 
and the estimate (|5.60l) are employed in the last one. 

Moreover, according to the solvability of L^-Dirichlet problems (see |32p Theorem 1.3]), it is not hard to 
see ||(Vp)*jjx, 2 (gQ) < C\\Vp\\]^ 2 (^Qf^y This together with (|5.74l) leads to 

II(V5)*IIl 2(9C2) < <^{ll5llL2(9n) + II-^II^(O)}- 
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Remark 5.8. Let Lq = —div(^V), and u G R”^) be a weak solution to Lq{u) = 0 in O. For any 

B{P, R) C AB C n, we may assume P = 0 and P = 1 from the translation and rescaling arguments. Then 
due to the interior regularity theory (see m Theorem 4.11]), we have ||u||Hfc(B) < C{k,m,d, fi)\\u\\]^ 2 (^ 2 B)i 
where P^(n;R™') = Moreover, it follows from the Sobolev imbedding theorem (for 2k > d) 

that |u(0)| < ||u||loo(_b) < C\\u\\jjk(^B) ^ C\\u\\i 2 ^ 2 B)- Note that v = also satisfy Lq{v) = 0 in fl. In 
fact, we thus have |V^u(0)| < C\\V‘^u\\i 2 (^ 2 B)- By Cacciopolli’s inequality (|2.14ji (or see [TSl Theorem 4.1]), 
we acquire |V^u(0)| < C\\Vu\\i 2 (^ 4 By Let u{x) = u{Rx) be the scaled function, a routine computation gives 

B VB(0,4iJ) ^ 

Remark 5.9. Let uq = v + p and u, p be given in (|5.56l) . As shown in the proof of Lemma 15.61 we can also 
prove the following results: 

l|V^no||L2(s^^) < ||V^u||i2(Rd) + ||V^/?||£2 (s^^) 

< C{\\F\\l2(^q) + ||5||L2(aQ)} +C'e“2||(Vp)*||i2(gt^) < Ce'z |||f||^2(q) + ||ff]|L2(9Q)}, 

(5.75) 

where A^ is a positive integer. Note that we use the Calderon-Zygmund theorem (see [181 Theorem 7.22]) 
for V and the estimate (15.671) in the second inequality, and Holder’s inequality in the last one. Also, 

ll^o||Hi(n\Sjve) - (5.76) 

which is the counterpart of the estimate (15.641) . Note that if we combine (15.271) . (15.751) and (15.761) by setting 
N = A, then we have 

||u£ - Uo - eXO,eS£{'lp4eUo) - £Xk,6Sei'll^4£'^kUo)\\< Cs^ {\\F\\l 2 (^q) + 1131^2(90)} (5.77) 

under the assumption that 11 is just a bounded Lipshitz domain. The estimate (15.771) is similar to (15.481) in 
Lemma 15.61 but the estimate (I5.48P is sharp in the sense of the integrability of the given data F, and it is 
also the reason why we mollify 'i/’ 2 e(Vuo) and ijj^eUo) twice in Lemma 15.61 

Theorem 5.10. Let Q be a bounded Lipschitz domain. Suppose that the coefficients of satisfy (HID-dLSl), 
and additional condition A = A*. Assume Ue G 77^(11;R™) is the weak solution to the Neumann problem 
= B in Q, and = g on dLl, where F G LP(H;R"*) forp = and g G L‘^{dLl;'SF). Then 

1 /* ^ 

|- / {\ue\^ + \Vue\‘^)dxy <C{\\FUpin) + \\9\\LHan)}, (5-78) 

£ 7n\Se 

where C depends only on p, K,m,d,p and H. 

Proof. Without loss of generality we may assume that l|Pl|LP(f 2 ) + || 5 ||L 2 (gf^) = 1. Let 

We = Ue - Uo - eXO,eSe{'^ieUo) “ £Xk,6Se{f^4e'^kUo) 

Hence, it follows from (15.481) and (|5.64p that 

lkellHi(r2\Se) — lko||//i(f^\Ss) + ll'“^£llHi(r2) < Ce^. 

Note that S‘^(^p 4 eUo) and S‘^{'ilj 4 eVkUo) is supported in therefore their 77^-norms vanish on H \ The 
proof is completed. □ 

Remark 5.11. If we additionally assume A G C'^(R'’*) in Theorem 15.101 and P = 0. Then, the estimate 
p5.78p in fact leads to (the Rellich estimate) ||Vue||^ 2 (gQ) < C\\g\\i 2 (^QQp where C is independent of e. The 
proof will be given in another place. Also, by referring to [291 Remark 3.1], the reader can prove it without 
real difficulties. 
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Proof of Theorem II.3L (i) With the help of the preceding Lemma 15.51 we can now prove the first part 
of Theorem 11.31 Note that 

\\Ue - Wo||l 2 (q) < \\Ue - Uq - eXO,eSe{lp46Uo) - eXj,eSs{i^4e'^jUo)\\L 2 (^^) 

+ £\\XO,sSei'll^4.eUo)\\L\n) + £:\\Xj,eSeilp4e'^ jUo)\\L2(^^^ 

< Cs{\\F\\l2^q) + ||5llBi/2,2(gn)} + C'e||uo||Hi(r2) 

< Cs{\\F\\l2^q-^ + \\g\\Bi/2,2(^gQ)}, 

where we use the estimates (I2.21h and (15.321) in the second inequality. In the last one, we mention that 

(ii) The second part of Theorem 11.31 follows from Lemma 15.61 As shown before, we have 

he - ^^o||l2(q) < \\Ue -Uq- £X0,eSl{'4^2eUQ) - £Xj,eSl{h2e^ jUo)\\L^{n) 

+ hXQ,eSl{llj2eUo)\\L^(n) + £\\Xj,eSl{llj2e'^ jUQ)\\L^(Q) 

< C'e 2 |||_F||^^^^^ + 1151^2(90)} + C'e{||5e(T/'2eUo)||L2(R'i) + ||S'e(V’ 2 £Vuo)||i;, 2 (]Rd)} 

< + WahHon) + + ||Vno||^^^^J 

<Ce5{||Fl|^^^^) + ||5|U2(ao)}. 

In the second inequality, we employ the estimates (12.211) and (15.481) . In the third one, we use the estimate 

( |2-23p . For the last one, we note that ||uo||wi.p(o) <Cho\\m(n) < C'{||-^IIlp(o) +1151^2(90)}, where p = ^. 

By duality argument, we will prove the following estimate 


he - < C£\\uo\\H2(^fi), 


(5.79) 


Let We be defined in ()5.35l) . For any <I> G L^+i (0; M™), let (/>£, 4>o be the solutions of (j5.33p (the existence 
is given by Lemma l3.Ij) . Then it follows from Lemma 15.61 that 


\\(j)e-4>0- £Xo,e^ehl2e(po) “ eXfc,£5'2(V'l2£Vfc4) < Ce2 ||$|| ^ 


L3+T(o) 


where Xk correctors associated with £* for fe = 0, • • • ,d. 

Proceeding as in the proof of Lemma 15.51 we arrive at (j5.45l) , which is 


/ 


We^dx 


< 


C'|||^^0||Hl(O\E9e)ll<^£|lfl'l(f^W9e) + ^ll'^O ||//2 (q) ||(()£ || jyl |. 


Due to Lemma 13.11 and Holder inequality, it is easy to see 


- " "L3+2(n) “ " "L3+T(n) 


(5.80) 


(5.81) 


(5.82) 


Set ^e = 4>e-4>0- £Xo,e^hh2e(po) “ £Xk,e^hh2e'^k4>o) , and then we have 


- ll'^£llHl(f2\S9e) + 11 <^0 || (t2\S9e) - ^ ^ ^ ^ ^ ^ i ^ (Q) 


(5.83) 


where we use the estimates (j5.76p and (j5.80l) in the second inequality. Collecting ()5.8ip . (j5.82p . (|5.83p . we 
obtain 


/ 

Jn 


We^dx 
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where we use the estimate (j5.29j) in the last inequality. This implies 

< Ce\\uo\\H2(n). 

and by the same procedure as we did in (i), it is not hard to derive the estimate (|5.79p . where we need to 
employ the estimate (|2.2ip for p = and ||Xfc||i,p(y) < C\\xk\\H^(Y) with k = 0, - ■ ■ ,d, (due to Sobolev 
embedding theorem). The details are left to readers, and we complete the proof. □ 
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